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TEMPORAL JOINT SEARCHER AND CHANNEL 

ESTIMATORS 

BACKGROUND 



[0001] This application is related to the following United States Patent applications, all 

5 simultaneously filed herewith: United States Patent Application Serial 10/ , 

(attorney docket 2380-776), entitled "Multi-Dimensional Joint Searcher And Channel 

Estimators"; United States Patent Application Serial 10/ , (attorney docket 2380- 

795), entitled "Spatial Joint Searcher And Channel Estimators"; United States Patent 

AppUcation Serial 10/ , (attomey docket 2380-797), entitled "Spatio-Temporal 

10 Joint Searcher And Channel Estimators", all of which are incorporated by reference 
herein. 

[0002] FIELD OF THE INVENTION 

[0003] The present invention pertains to wireless telecommunications, and particularly 
to apparatus and method for determining a channel estimate for use in reconstructing 
15 data symbols transmitted over a channel. 

[0004] RELATED ART AND OTHER CONSIDERATIONS 

[0005] A wireless telecommunications unit typically includes both a transmitter and 
receiver for communicating with other wireless telecommunications units over a 
communication link. For wireless communications, the communication link is typically 

20 over an air interface (e.g., radio frequency interface). As used herein, a "wireless 
telecommunications unit" with its "wireless telecommunications receiver" can be 
included in a network node (e.g., a radio access network node such as a base station 
node, also called Node-B) or a terminal. Such "terminals" include mobile terminals 
such as user equipment units (UEs), which have also been called mobile stations, and 

25 include by way of example mobile telephones ("cellular" telephones), laptops with 
mobile termination. Thus, terminals can be, for example, portable, pocket, hand-held. 
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computer-included, or car-mounted mobile devices which communicate voice and/or 
data with radio access network. Altematively, the terminals can be fixed wireless 
devices, e.g., fixed cellular devices/terminals which are part of a wireless local loop or 
the like. 

5 [0006] As shown simply in Fig. 32, a wireless teleconraiunications system includes a 
transmitting antenna 2300T and a receiving antenna 2300R. Channel 2302 describes 
the relation between the transmitting antenna 2300T and the receiving antenna 2300R, 
including the wireless interface. A signal, typically modulated into pulses, is 
transmitted over channel 2302 from transmitting antenna 2300T to receiving antenna 

10 2300R. The signal can comprise a "s3niibol" or a string of series of symbols, depicted 
as "m" in Fig. 32. The signal can carry user data and/or certain control data (e.g., a 
pilot bit or pilot sequence). The signal m as transmitted by the transmitting antenna 
2300T is convoluted with a channel impulse response h of the channel, so that the 
received signal at the receiving antenna 2300R is m * /i (e.g., m convoluted with h). 

15 The received signal m his applied to base band processing functionality 2304 of the 
receiver where the received signal undergoes radio frequency processing. The data 
portions of the received signal are applied to a detector 2306, which may be, for 
example, a demodulator such as a RAKE receiver. 

[0007] Most modem detectors attempt to recover a symbol estimate m from the 
20 received signal m^h . To do so, most sophisticated detectors expect to receive a 
"channel estimate" for use in modeling the channel over which the signal was 
transmitted. The accuracy of this channel estimate influences the accuracy and 
performance of the detector in estimating the actual symbol received over the channel. 

[0008] The modeling of the channel (which is necessary for most detectors) is 
25 facilitated by the control data, often in the form of a pilot bit or pilot sequence, which is 
transmitted by the transmitter. The control data, hereinafter referenced as "pilot data" 
for simplicity, is of a known or recognizable format or pattern. The pilot data is 
typically transmitted periodically by the transmitter source, and thus receipt of 
repetitions of the pilot data can be expected at the receiver at successive intervals. In 
30 view of factors such as relative motion of the transmitter and receiver, the successive 
intervals are not necessarily constant. The pilot data can be transmitted simultaneously 
with, or otherwise interspersed with, the user data. 



3 

[0009] In order to utilize the pilot data, wireless receivers typically include both a 
searcher and a channel estimator, such as searcher 2308 and channel estimator 2310 
shown in Fig. 32. For control data, the received signal m ^ his applied to the searcher 
2308, which determines a time of arrival (TOA). The time of arrival is then applied to 
5 the channel estimator 2310, which uses the time of arrival to determine the channel 
estimate h and then provides the channel estimate h to the detector 2306. Using the 
channel estimate h , the detector develops its estimate of the symbol, e.g., m . 

[00010] The receiver may receive the an original signal (e.g., short pulse signal) from 
the transmitter source through open space over a single, direct propagation path. 
10 Alternatively, in another environment having obstacles or other surfaces, the receiver 
may receive the same original signal over multiple propagation paths. In the multiple 
path case, the received signal appears at the receiver as a stream of pulses, each pulse 
having a different time delay in view of the corresponding propagation multipath over 
which the signal travelled, as well as possibly different amplitude and phase. 

15 [0001 1] Multipaths are created in a mobile radio channel by reflection of the signal 

from obstacles in the environment such as buildings, trees, cars, people, etc. Moreover, 
the mobile radio channel is dynamic in the sense that it is time varying because of 
relative motion affecting structures that create the multipaths, or due to movement of 
structures and objects in the surroundings (even if the transmitter and receiver are 

20 fixed). For a signal transmitted over a time varying multipath channel, the received 
corresponding multiple paths vary in time, location, attenuation, and phase. 

[00012] Some wireless telecommunications receivers capitalize upon the existence of 
the multipaths in order to achieve various advantages. Such receivers typically operate 
on the baseband signal to search for and identify the strongest multipaths along with 

25 their corresponding time delays. The receiver has a filter which operates on a power 
delay profile of the signal. The power delay profile can be conceptualized as a time- 
averaged refinement or other derivation of the channel impulse response. The searcher 
attempts to locate peaks in the power delay profile, each peak corresponding to arrival 
of a wavefront of the signal from a respective multipath. In many searchers the peeiks 

30 also correspond to a channel tap of the filter. 
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[00013] A channel estimate h as applied to the detector therefore comprises a set of 
both time of arrivals (TO A) and complex channel coefficients, each pair of TOA and 
channel coefficients being associated with one of the arriving wavefronts. In other 
words, each arriving wavefront has a pair of members in the set, e.g., a TOA and a 
5 channel coefficient. The channel coefficients thus actually form a channel impulse 
response vector, so that the terms "channel coefficient" and "channel coefficients" as 
used hereinafter should be understood to refer to a channel impulse response vector. If 
there is only one wavefront, there is only one TOA and one channel coefficient in the 
set (one channel coefficient in the channel impulse response vector) . But for plural 
10 arriving wavefronts, there are a corresponding plurality of TOAs and channel 

coefficients. Ideally, the channel estimate h should provide as good an estimate of the 
channel impulse response as possible, thereby increasing performance of the detector as 
the detector makes its estimate m of the transmitted symbol m. 

[00014] The channel estimate is then supplied to the detector, such as RAKE type of 
15 demodulators. A RAKE demodulator typically allocates a number of parallel 

demodulators (called RAKE fingers) to the strongest multipath components of the 
received multipath signal as determined by the multipath search processor. In a 
wideband code division multiple access (WCDMA) radio access network, the outputs 
of each of the RAKE fingers are diversity-combined after corresponding delay 
20 compensation to generate a "best" demodulated signal that considerably improves the 
quality and reliability of the radio communications system. 

[00015] Conventionally, wireless telecommunications receivers first use their searchers 
to ascertain time of arrival of a wavefront. Subsequently, after the time of arrival has 
been determined by the searcher, the channel estimator utilizes the time of arrival to 
25 calculate a channel coefficient, which expresses both amplitude and phase of the signal. 

[00016] Some wireless telecommunications units have more than one antenna for 
receiving a same signal. In the prior art, the searcher attempts to locate peaks in the 
power delay profile for each antenna separately. In other words, for each antenna the 
searcher works more or less independently. See, for example, US Patent Publication 
30 US 2002/0048306, which is incorporated herein by reference. As such, the prior art 
searchers are essentially one dimensional. 
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[00017] As indicated above, the performance of a wireless receiver is considerably 
dependent upon the accuracy of the peak determination, i.e., time of arrival 
determination, performed by the searcher. The better the peak determination of the 
searcher, the better will be the overall performance of the receiver (e.g., less error rate). 
5 But in many instances it may be difficult for a searcher to find an actual peak in a 
power delay profile. As mentioned previously, in m£my searcher algorithms the peak 
corresponds to a channel tap. With such difficulty there is considerable risk of 
incorrectly choosing a peak. Moreover, it can then be difficult to estimate the actual 
channel tap value. Channels with low signal to noise ratios (SINRS) are particularly 
10 susceptible to these difficulties. 

[00018] What is needed, therefore, and an object of the present invention, is provision 
of apparatus and method for providing an improved channel estimate for a wireless 
telecommunications receiver. 

BRIEF SUIMIMARY 

[00019] A wireless conmiunication receiver comprises a joint searcher and channel 
estimator which provides a channel estimate which can take into consideration a 
frequency shift (e.g., doppler shift) occasioned, e.g., by relative movement of a 
transmitter and the receiver or movement of a signal path-affecting object. In providing 
the channel estimate, the joint searcher and channel estimator essentially concurrently 
considers plural signals received by an antenna element of the array, the plural signals 
being obtained from a series of successive sets of pilot data as detected by the antenna 
element. The time of arrival and the channel coefficient are essentially concurrently 
determined by the joint searcher and channel estimator. The joint searcher and channel 
estimator applies the channel coefficient and the time of arrival to a detector which 
provides, e.g., a symbol estimate. 

[00020] The wireless communication receiver can be either a mobile terminal or a 
network node (e.g., a radio access network node such as a base station node, also called 
Node-B). The relative movement of the transmitter and receiver which gives rise to the 
doppler shift (which is factored into the determination of the channel estimate by the 
30 joint searcher and channel estimator) can occur by reason, e.g., of movement of a 

mobile terminal. Such movement can be presumed to be at a relatively constant speed. 



20 
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particularly with reference to the time base over which the series of successive sets of 
pilot data are processed by the joint searcher and channel estimator. Moreover, the 
filtering can be done adaptively, e.g., if the mobile station does not move, a longer 
period of data can be utilized than if the mobile station does not move. 

5 [00021] For a single antenna element, the signals over respective sampling windows for 
each of the series of successive sets of pilot data are contemporaneously processed in 
order to determine both time of arrival and the channel coefficient. Accordingly, the 
joint searcher and channel estimator is considered a two dimensional unit, A first 
dimension is with reference to a time index of the sampling window employed for each 
10 of the sets of pilot data, i.e., a sampling window time index. 

[00022] A second dimension is a temporal dimension imparted by the time interval 
reflected by the successive sets of pilot data. This temporal dimension, which involves 
essentially simultaneous and concurrent processing together of signals for each of the 
successive sets of pilot data in order to determine the time of arrival and channel 
15 coefficient, bestows on the joint searcher and channel estimator the distinction of being 
a "temporal" joint searcher and channel estimator. 

[00023] The temporal joint searcher and channel estimator may take differing 
embodiments and have differing implementations. In one example, illustrative 
embodiment the temporal joint searcher and channel estimator includes a non- 
20 parametric type correlator (e.g., a correlator which performs a Fast Fourier Transform 
(FFT) calculation). In another example, illustrative embodiment the temporal joint 
searcher and channel estimator utiUzes a parametric approach. 

[00024] Concurrently using signals from each of the successive sets of pilot data, the 
joint searcher and channel estimator looks for pilot data in a sampling window for each 
25 of the successive sets of pilot data, and thereby factors a possible doppler shift into the 
generation of the time of arrival and channel coefficient for each wavefront having a 
peak as seen in the sampling windows. 

[00025] In so doing, for each sampling window the joint searcher and channel estimator 
stores the convoluted baseband signals obtained from the antenna. One example way of 
30 storing the signals for a sampling window is in a matrix, e.g., an antenna signal matrix. 
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In constructing the antenna signal matrix, each of the sets of pilot data to be considered 
together by the joint searcher and channel estimator is represented by a doppler 
frequency index, or "doppler index". The joint searcher and channel estimator stores in 
the antenna signal matrix a complex value indicative of the signal received in the 
5 sampling window for each set of pilot data. The position or location of the complex 
value indicative of the signal received is determined by two indexes. The first index, 
conceptualized as being along the X axis of the antenna signal matrix, is the sampling 
window time index. The sampling window time index points to a time in the respective 
sampling window relative to a start of the sampling window. The second index, 
10 conceptualized as being along the Y axis of the antenna signal matrix, is the doppler 
frequency index. 

[00026] In the embodiment in which the joint searcher and channel estimator includes a 
correlator which performs a Fast Fourier Transform (FFT) calculation, the correlator 
considers a dimensional receptivity vector. The dimensional receptivity vector includes 

15 values of a comparably numbered sampling window time index for differing sets of 

pilot data. In other words, the dimensional receptivity vector is a column of the antenna 
signal matrix. A rotation speed, or frequency, of the dimensional receptivity vector for 
the sampling window time instance, indicated by differing imaginary parts (9 values) 
for the members of the dimensional receptivity vector, reflects a moving or rotation of 

20 the complex baseband caused, e.g., by the movement of the transmitter and receiver 

(e.g., of the mobile terminal) or movement of a signal-path affecting object or structure. 
Each dimensional receptivity vector represents one of plural possible doppler 
shifts/frequencies. In conjunction with the Fast Fourier Transform (FFT) calculation, 
the correlator calculates the doppler frequency as Y(n,t) = FFT(n,X(:,t)), wherein t is 

25 the sampling window time index; X(n,t) is the complex antenna matrix; and, n is the 
doppler frequency index. For a CDMA receiver, the correlator calculates Y(n,t) = S Cj 
*FFT(n,X(:,t)), j = 1,K, wherein Cj is a coding sequence symbol value j and K is a 
length of the coding sequence. 

[00027] In the Fast Fourier Transform (FFT) calculation, the correlator output 
30 comprises Y(n,t). An analyzer of the joint searcher and channel estimator determines a 
maximum absolute value |Y(n,t)|niax from the correlator output Y(n,t). A sampling 
window time index t_max at which | Y(n,t)|max occurs is chosen as the time of arrival of 
the arriving wavefront; a frequency index n_max at which |Y(n,t)|niax occurs is chosen 
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as the doppler shift frequency of the arriving wavefront. An amplitude for the arriving 
wavefront is chosen by dividing |Y(n,t)|max by the number of pilot data sets considered. 

[00028] In another embodiment the joint searcher and channel estimator includes a 
parametric estimator which generates a parametric estimation output vector which is 
5 utilized by a channel estimate generator to generate the time of arrival and the channel 
coefficient. The parametric estimation output vector includes temporal frequency value 
and temporal amplitude value for each time instant. The parametric estimation output 
vector has a sampling window time index and temporal frequency parameter value for 
each time index. The channel estimate generator uses spatial amplitude values of 

10 elements of the parametric estimation output vector to determine the time of arrival and 
the doppler shift frequency of the arriving wavefront. Wavefronts in the sampling 
windows are associated with each element of the parametric estimation output vector 
which has the sufficiently high absolute value. The channel estimate generator uses a 
sampling window time index for an element of the parametric estimation output vector 

15 having a sufficiently high absolute value as the time of arrival of the corresponding 
arriving wavefront. The doppler shift of the arriving wavefront is the temporal 
frequency parameter value of the identified time of arrival. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[00029] The foregoing and other objects, features, and advantages of the invention will 
20 be apparent from the following more particular description of preferred embodiments as 
illustrated in the accompanying drawings in which reference characters refer to the 
same parts throughout the various views. The drawings are not necessarily to scale, 
emphasis instead being placed upon illustrating the principles of the invention. 

[00030] Fig. 1 is a schematic view of an example, generic wireless telecommunications 
25 receiver which includes a joint searcher and channel estimator. 

[00031] Fig. 2 A and Fig. 2B are schematic views of differing example embodiments of 
spatial joint searcher and channel estimators, each shown with an antenna array. 
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[00032] Fig. 3 is a diagrammatic view illustrating a signal emanating from a 
transmitting antenna along three separate multipaths to an antenna array of a wireless 
telecommunications receiver. 

[00033] Fig. 4 is a diagranraiatic view of a wavefront travelling toward an antenna 
5 array. 

[00034] Fig. 5 A and Fig. 5B are diagrammatic views depicting signals obtained upon 
arrival of a wavefront at an antenna array. 

[00035] Fig. 6 is a diagrammatic view of an antenna signal matrix. 

[00036] Fig. 7 is a flowcharting showing representative basic steps performed by a 
10 matrix analyzer and channel estimate generator of an example embodiment of a spatial 
joint searcher and channel estimator, with the matrix analyzer using a non-parametric 
analysis technique. 

[00037] Fig. 8A, Fig. 8B, Fig. 8C(1), Fig. 8C(2), and Fig. 8C(3) and are diagrammatic 
views depicting results of a comparative operational evaluation contrasting 
15 performance of a spatial joint searcher and channel estimator with a conventional 
searcher. 

[00038] Fig. 9A is a diagrammatic view of an antenna signal matrix; an antenna weight 
vector; and a non-parametric output estimation vector. 

[00039] Fig. 9B is a diagrammatic view of an antenna signal matrix and a parametric 
20 output estimation vector. 

[00040] Fig. 10 is a flowcharting showing representative basic steps performed by 
matrix analyzer and channel estimate generator of an example embodiment of a spatial 
joint searcher and channel estimator, with the matrix analyzer using a parametric 
analysis technique. 

25 [00041] Fig. 1 1 is a diagrannmatic view illustrating a coherent combination of signal 
outputs by a joint searcher and channel estimator. 
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[00042] Fig. 12 is a schematic view for illustrating how an antenna weight vector 
facilitates the coherent combination illustrated in Fig. 11. 

[00043] Fig. 13A is a schematic view of an example embodiment of a temporal joint 
searcher and channel estimator shown with an antenna array, the temporal joint 
5 searcher and channel estimator comprising a matrix analyzer which employs a non- 
parametric analysis technique. 

[00044] Fig. 13B is a schematic view of an example embodiment of a temporal joint 
searcher and channel estimator shown with an antenna array, the temporal joint 
searcher and channel estimator comprising a matrix analyzer which employs a 
10 parametric analysis technique. 

[00045] Fig. 14 is a diagrammatic view depicting a sequence of sets of pilot data and 
user data received by a receiver which utilizes a temporal joint searcher and channel 
estimator, as well as an antenna signal matrix utilized by the temporal joint searcher 
and channel estimator. 

15 [00046] Fig. 15 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a 
temporal joint searcher and channel estimator, with the matrix analyzer using a non- 
parametric analysis technique. 

[00047] Fig. 16A is a diagrammatic view of an antenna signal matrix; a doppler weight 
20 vector; and a non-parametric estimation output vector for a temporal joint searcher and 
channel estimator. 

[00048] Fig. 16B is a diagrammatic view of an antenna signal matrix and a parametric 
estimation output vector for a temporal joint searcher and channel estimator. 

[00049] Fig. 17 is a flowcharting showing representative basic steps performed by a 
25 matrix sinalyzer and channel estimate generator of an example embodiment of a 
temporal joint searcher and channel estimator, with the matrix analyzer using a 
parametric analysis technique. 
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[00050] Fig. 18A is a schematic view of an example embodiment of a spatio-temporal 
joint searcher and channel estimator shown with an antenna array, the spatio-temporal 
joint searcher and channel estimator comprising a matrix analyzer which employs a 
non-parametric analysis technique. 

5 [00051] Fig, 18B is a schematic view of an example embodiment of a spatio-temporal 
joint searcher and channel estimator shown with an antenna array, the spatio-temporal 
joint searcher and channel estimator comprising a matrix analyzer which employs a 
parametric analysis technique. 

[00052] Fig. 19 is a diagrammatic view depicting a sequence of sets of pilot data and 
10 user data received by a receiver which utilizes a combined spatial/temporal joint 
searcher and channel estimator, as well as an antenna signal matrix utilized thereby. 

[00053] Fig. 20 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a spatio- 
temporal joint searcher and channel estimator, with the matrix analyzer using a non- 
15 parametric analysis technique. 

[00054] Fig. 21 is a diagrammatic view of an antenna signal matrix; a doppler weight 
and antenna weight vector; and a non-parametric estimation output vector for an 
example embodiment of a spatio-temporal joint searcher and channel estimator which 
operates in a three dimensional essentially concurrent mode. 

20 [00055] Fig. 22A and Fig. 22B are diagrammatic views depicting operation of a first 
altemative implementation of a non-parametric, sequential spatio-temporal joint 
searcher and channel estimator. 

[00056] Fig. 23 describes the procedure of non-parmetric approach for spatio-temporal 
sequenced method where the spatial processing is followed by the temporal processing. 

25 [00057] Fig. 24 A and Fig, 24 Bare diagrammatic views depicting operation of a second 
altemative implementation of a non-parametric, sequential spatio-temporal joint 
searcher and channel estimator. 
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[00058] Fig. 25 describes the procedure of non-parmetric approach for spatio-temporal 
sequenced method where the temporal processing is followed by the spatial processing. 

[00059] Fig. 26 is a diagrammatic view of an antenna signal matrix and a parametric 
estimation output vector for an example embodiment of a spatio-temporal joint searcher 
5 and channel estimator. 

[00060] Fig. 27 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a spatio- 
temporal joint searcher and channel estimator, with the matrix analyzer using a 
parametric analysis technique. 

10 [00061] Fig. 28 A and Fig. 28B are diagranunatic views depicting operation of a first 
altemative implementation of a parametric, sequential spatio-temporal joint searcher 
and channel estimator. 

[00062] Fig. 29 describes the procedure of the parametric approach for spatio-temporal 
sequenced method where the spatial processing is followed by the temporal processing. 

15 [00063] Fig. 30A and Fig. BOB are diagrammatic views depicting operation of a second 
altemative implementation of a parametric, sequential spatio-temporal joint searcher 
and channel estimator. 

[00064] Fig. 31 describes the procedure of the parametric approach for spatio-temporal 
sequenced method where the temporal processing is followed by the spatial processing. 

20 [00065] Fig. 32 is a schematic view of a conventional wireless telecommunications 
receiver. 



DETAILED DESCRIPTION OF THE DRAWINGS 

[00066] In the following description, for purposes of explanation and not limitation, 
specific details are set forth such as particular architectures, interfaces, techniques, etc, 
25 in order to provide a thorough understanding of the present invention. However, it will 
be apparent to those skilled in the art that the present invention may be practiced in 
other embodiments that depart from these specific details. In other instances, detailed 
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descriptions of well-known devices, circuits, and methods are omitted so as not to 
obscure the description of the present invention with unnecessary detail. Moreover, 
individual function blocks are shown in some of the figures. 

[00067] Fig. 1 shows an example, generic wireless telecommunications receiver 20 
5 which, as mentioned before, can be included in a network node or a terminal, e.g. 
mobile terminal. The wireless telecommunications receiver 20 includes an antenna 
structure or array 22; a joint searcher and channel estimator 24; a detector 26; and, a 
timing and control unit 28. Optionally, as depicted by broken line, the receiver 20 may 
include a code sequence generator 30. 

10 [00068] As broadly employed herein, the antenna array 22 can comprise one or more 
antenna elements. Signal(s) from the antenna array 22 are applied both to joint searcher 
and channel estimator 24 and detector 26. The signal(s) from the antenna array 22 
comprise a channel impulse response vector if the antenna array 22 comprises more 
than one antenna element. 



15 [00069] In the likely event that the signal(s) have been encoded by, e.g., a spreading 
code or the like, both joint searcher and channel estimator 24 and detector 26 are 
connected to operate in conjunction with code sequence generator 30. The timing and 
control unit 28 generates timing (e.g., synchronization) and control signals which are 
provided to detector 26 and to joint searcher and channel estimator 24. 



20 [00070] It will be appreciated that the receiver may include, e.g., downstream from the 
antenna array, certain radio frequency processing functionality and radio frequency 
demodulating functionality, so that the signals applied to joint searcher and channel 
estimator 24 and detector 26 are baseband signals. The illustrated structure of wireless 
telecommunications receiver 20 of Fig. 1 thus essentially concerns processing of the 

25 baseband signal(s). 

[00071] Various non-limiting, representative examples of differing embodiments of 
joint searcher and channel estimators are described below. Ensuing descriptions of 
operation of wireless telecommunications receivers with these differing embodiments 
are premised on certain assumptions. Some of these assumptions are related to a 
30 channel model which conceptualizes electromagnetic fields as arriving in a discrete 
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number of wavefronts at the wireless telecommunications receiver, and particularly 
arriving at one or more antenna elements which may be employed in the antenna array 
22. 

[00072] As used herein, a "sampling window" comprises consecutive time slots (or, in 
5 a CDMA system, for example, "chips") obtained from a given antenna and analyzed by 
a joint searcher and channel estimator. As described in more detail hereinafter, 
embodiments of joint searcher and channel estimators operate upon an antenna signal 
matrix formed from plural sampling windows. In some embodiments, hereinafter 
referenced as "spatial" joint searcher and channel estimators, the antenna signal matrix 

10 is formed from sampling windows obtained from plural antennas. In other 
embodiments, hereinafter referenced as "temporal" joint searcher and channel 
estiniators, the antenna signal matrix is formed with respect to a single antenna, but 
formed from sampling windows obtained by that antenna for successive sets of pilot 
data (occurring over time). In yet other embodiments, hereinafter referenced as spatio- 

15 temporal joint searcher and channel estimators, the antenna signal matrix is formed 
both spatially and temporally. 

[00073] For purposes of the technology described herein, the antenna array 22 is 
conceptualized as acquiring "dimensionally differentiated" signals. The joint searcher 
and channel estimator essentially concurrently uses the dimensionally differentiated 

20 signals provided by the antenna array for determining, for each arriving wavefront, both 
a time of arrival (TO A) and a channel coefficient. For the spatial joint searcher and 
channel estimator, wherein the antenna structure comprises an array of plural antennas 
having spaced apart or spatially separated antenna elements, the signals acquired by 
different antennas of the array are dimensionally differentiated with regard to a spatial 

25 dimension. For the temporal joint searcher and channel estimator, wherein the antenna 
structure comprises an antenna which provides signals for each of successive sets of 
pilot data received at separated time intervals, the signals acquired by the antenna are 
dimensionally differentiated with regard to a temporal or time dimension. For the 
spatio-temporal joint searcher and channel estimator, having both the antenna structure 

30 comprising an array of plural antennas and one or more antennas receiving the 

successive sets of pilot data, the signals acquired by the antenna are dimensionally 
differentiated with regard both to a spatial dimension and a temporal or time dimension. 
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[00074] The joint searcher and channel estimators are, in some instances, said to 
perform a "concurrent" determination of time of arrival £ind some other quantity, e.g., 
direction of arrival or doppler shift frequency. In this sense "concurrent" means that the 
quantities or determinations could be derived in parallel from a result of an outcome- 
5 determinative operation, e.g. a non-parametric technique such as a Fast Fourier 
Transform or a parametric technique. 

[00075] SPATIAL JOINT SEARCHERS/ESTIMATORS 

[00076] In some embodiments, the joint searcher and channel estimator 
contemporaneously processes the signals from plural antennas over a sampling window 

10 in order to determine both time of arrival and the channel coefficient. In these 

embodiments the joint searcher and channel estimator is essentially a two dimensional 
unit, with a second dimension being a spatial dimension imparted by the spacing of the 
plural antennas of the array. This spatial dimension, which involves essentially 
simultaneous and concurrent processing of signals from the plural antennas for the 

15 array in order to determine the time of arrival and channel coefficient, bestows on these 
embodiments of the joint searcher and channel estimator the distinction of being a 
"spatial" joint searcher and channel estimator. 

[00077] The spatial joint searcher and channel estimator may take differing 
embodiments and have differing implementations. In one example, illustrative 
20 embodiment the joint searcher and channel estimator includes a non-parametric type 
correlator (e.g., a correlator which performs a Fast Fourier Transform (FFT) 
calculation). In another example, illustrative embodiment the joint searcher and 
channel estimator utilizes a parametric approach 

[00078] Fig. 2A illustrates one example embodiment of a spatial joint searcher and 
25 channel estimator 24-2A which uses a non-parametric technique for determining time 
of arrival and channel estimate, as well as an associated example antenna array 22-2A. 
The antenna array 22-2A includes, by way of non-limiting example, four antenna 
elements 22-2A-1 through 22-2A-4. While the antenna elements 22-2A-1 through 22- 
2A-4 are shown as forming a uniform linear array (ULA), it should be understood that 
30 antenna configurations other than a uniform linear are possible, and that the number of 
antenna elements in the antenna array may vary (e.g., the number of antenna elements 
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is not limited to four). 

[00079] There are coherency requirements for the antenna elements of antenna array 
22-2A, and for antenna elements for all other plural antenna arrays described herein. 
The coherency requirement can be fulfilled by the plural antenna elements being 
5 synchronized. Alternatively, even if the plural antenna elements are not synchronized, 
but their phase differences are known, the coherency requirement can be fulfilled by 
compensating for the known phase difference. 

[00080] The complex baseband signals obtained from the antenna elements are each 
applied to joint searcher and channel estimator 24-2 A, as well as to a detector (not 

10 illustrated in Fig. 2A). The joint searcher and channel estimator 24-2A comprises an 
antenna signal matrix handling unit 40-2A. In one particular example manifestation, 
antenna signal matrix handling unit 40-2A includes antenna signal matrix generator 42- 
2A and antenna signal matrix memory 44-2A. A matrix analyzer, which for the non- 
parametric technique of Fig. 2 A can be a correlator 50-2 A, operates on complex values 

15 stored in anterma signal matrix memory 44-2A. The correlator 50-2A preferably 

comprises a filter. The correlator 50-2A generates certain output values, which may be 
stored, e.g., in correlator output value memory 52-2 A. The joint searcher and channel 
estimator 24-2 A further comprises a channel estimate (CE) generator 60-2 A. In the 
illustrated example embodiment, the channel estimate (CE) generator 60-2A comprises 

20 a correlator output analyzer 62-2A and a detector interface 64-2A. The detector 
interface 64-2A generates, for each wavefront, both a time of arrival (TOA) and a 
channel coefficient (CC). In Fig. 2A, the time of arrival and channel coefficient output 
by detector interface 64 are applied to the detector on lines 66-2A and 68-2A, 
respectively. 

25 [00081] In Fig. 2A, and other embodiments described herein, the transmitted 

electromagnetic signal is assumed to arrive at the receiver in a number of discrete 
electromagnetic wavefronts. A number of discrete electromagnetic wavefronts is 
presumed in order to accommodate the multipath phenomena discussed above. For 
example. Fig. 3 illustrates a signal emanating from a transmitting antenna 70 along 

30 three separate multipaths Pi, P2, and P3 to antenna array 22. Each multipath has its 
individual amplitude, and accordingly has an associated complex number "a" of the 
baseband signal and a time delay x. For example, multipath Pi has associated complex 
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number ai and associated time delay ti; multipath P2 has associated complex number a2 
and associated time delay T2; and so forth. As illustrated in Fig. 3, multipath Pi is a 
relatively direct path between transmitting antenna 70 and antenna array 22; while 
multipath P2 and multipath P3 are reflected off obstacles 722 and 723, respectively. 
5 Thus, the time delay Xi for multipath Pi is shorter than the time delay T2 for multipath 
P2, which in turn is shorter than the time delay T3 for multipath P3. Similarly, barring 
other phenomena, it would be expected that the complex number ai for multipath Pi is 
greater than the complex number a2 for multipath P2, and so forth. 

[00082] For sake of discussion, the electromagnetic wavefronts are assumed to be plane 
10 ("planar") electromagnetic wave fronts, such as the single wavefront 76 illustrated in 
Fig. 4 as traveling toward the antenna array. In all embodiments described herein it 
should be understood that the wavefronts need not be planar wavefronts, but that any 
other known form of wavefront may be considered in similar manner. Moreover, it 
should be kept in mind that Fig. 4 represents arrival of only one wavefront, but 
15 typically plural wavefronts are incident on an antenna array. 

[00083] As further shown by Fig. 4, due to incidence of an individual wavefront, the 
output (e.g., signal) from each antenna element has its version of the complex number 
for the wavefront. For example, for a wavefront for the first multipath Pi of Fig. 3, the 
antenna element 22-1 outputs a complex number ai.i, antenna element 22-2 outputs a 

20 complex number ai.2, and so forth. The numbers are complex, and in the particular case 
that (1) the antenna elements are identical; (2) there is coherency, and (3) the plane 
wave has constant amplitude within the width of array, the absolute values of the 
numbers are the same. Furthermore, with respect to the same arriving wavefront, each 
antenna detects the arriving signal as having a phase. For example, for the wavefront 

25 for the first multipath Pi of Fig. 3, the output of antenna element 22-1 has a phase 9 1.1, 
the output of antenna element 22-2 has a phase 0i.2, and so forth. 

[00084] Signals obtained upon arrival of a wavefront at a uniform linear array (ULA) 
antenna array are illustrated both in Fig. 5A and in Fig. 5B. Fig. 5A particularly shows, 
for each of four antennas 22-1 through 22-4, plane wave propagation over the antenna 
30 elements for a fixed time (chip) index, and resulting respective output pulses 78 (e.g., 
output pulses 78i through 784). For each corresponding antenna, Fig. 5B shows the 
pulse as a complex number and with the argument,©, of the complex number. The 
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argument (9) corresponds to the phase of the received signal. The rate at which the 0 
values change (e.g., the rate at which the phase rotates) over time is known as the phase 
rotation speed, or frequency. The phase rotation for the wavefront with this array of 
antennas is depicted by the increasing angular value of 0 through the range of 0i, 02, 03, 
5 04, and thus the frequency is the rate of change of this angular value over time. The 
phase rotation speed is constant. The speed of the linear phase propagation is 
dependent on the direction of arrival (DOA) of the incident wavefront. 

[00085] In the joint searcher and channel estimator 24-2A of Fig. 2A, the antenna 
matrix handling unit 40-2A samples the complex baseband signals from each antenna 
10 element. Using the sampled complex baseband signals, antenna signal matrix generator 
42-2A generates an antenna signal matrix such as antenna signal matrix 80 illustrated in 
Fig. 6. The antenna signal matrix 80 may be stored in any convenient fashion, such as 
antenna matrix memory 44-2 A. 

[00086] The antenna signal matrix 80 is a two dimensional functionally dependent 
15 matrix. In other words, complex samples are stored in antenna signal matrix 80 as a 
function of two different indexes. For the antenna signal matrix 80 shown in Fig. 6, a 
first index is a sampling window time index, illustrated along the X axis of Fig. 6. For 
embodiments which utilize spreading codes or similar codes, the first index may be, for 
example, a chip index. Thus, the sampling window time index points to a time in the 
20 sampling window relative to a start of the sampling window. In the antenna signal 

matrix 80 of Fig. 6, a second index, shown along the Y axis, is an antenna index (which 
serves as a dimensional differentiation index). The antenna index points to a different 
row of the antenna signal matrix 80, each row being associated with a different antenna 
element in antenna array 22. Fig. 6 shows four rows in antenna signal matrix 80 for 
25 consistency with the previous examples of an antenna array comprising four antenna 

elements. It is reiterated, however, that the number of antennas in an antenna array, and 
thus the number of rows in antenna signal matrix 80 and the maximum value of the 
antenna index, can vary from receiver to receiver, and that the choice of four antenna is 
only illustrative for sake of example. 

30 [00087] The antenna signal matrix 80 is conceptualized as storing '^dimensionally 
differentiated" signals acquired from the antenna array. For the spatial joint searcher 
and channel estimator, wherein the antenna structure comprises an array of plural 



antennas having spaced apart or spatially separated antenna elements, the signals 
acquired by different antennas of the array are dimensionally differentiated with regard 
to a spatial dimension. That is, for a given column of antenna signal matrix 80, the 
values in each row are dimensionally differentiated in the sense that they are, acquired 
5 from different antenna elements which are separated in a spatial dimension in view of 
the separate physical placement of each antenna element with respect to other antenna 
elements of the array. 

[00088] For sake of simplicity, the complex values stored in antenna signal matrix 80, 
including the complex values obtained from the antennas, are not illustrated in Fig. 6. 
10 Such complex values would be illustrated in a third dimension, e.g., out of the plane of 
Fig. 6. The antenna signal matrix 80 includes both complex white noise and (for the 
sake of the present illustration) a complex sample for at least one wavefront (planar or 
other known shape). As stored in antenna signal matrix 80, the wavefronts have known 
phase (temporal, non-coherent detection), and are modulated code sequences. 

15 [00089] In conjunction with the antenna signal matrix 80 of Fig. 6, and particularly a 
WCDMA case in which spacing of the antenna elements in the antenna array is not too 
far apart, the plane wavefront arriving at the antenna array can be considered to arrive 
in the same sampling window time index (or chip index). 

[00090] The complex values stored for each column of the antenna signal matrix 80 of 
20 Fig. 6 can be conceptualized as a dimensional receptivity vector. That is, a dimensional 
receptivity vector is formed with respect to a single sampling window time instance and 
with complex values from each of the plural antennas of the antenna array. Each 
element taken from a unique row of antenna signal matrix 80 has a different phase in 
the manner of the differing 0 values illustrated in Fig. 5. As received by the differing 
25 antenna elements, for the spatial joint searcher and channel estimator the change in 
phase over time is the frequency for the dimensional receptivity vector. If the wave 
arrives e.g. straight ahead, the angles could be the same. The phase rotation speed, or 
frequency, of the dimensional receptivity vector, for the sampling window time 
instance can be interpreted as the direction of arrival (DOA). Thus, each dimensional 
30 receptivity vector corresponds to a separate direction of arrival. There are plural 

possible frequencies for the dimensional receptivity vector, each of the plural possible 
frequencies corresponding to a different possible direction of arrival (DOA) of a 



20 



wavefront. For the non-parametric techniques herein employed, the plural possible 
frequencies can be a continuous range of frequencies. For sake of differentiating the 
plural possible frequencies, the plural possible frequencies are each represented by a 
frequency index. 

5 [00091] The channel estimate generator 60-2 A (see Fig, 2A) seeks to develop a 
"composite" channel estimate based on the complex values stored in antenna signal 
matrix 80. At this point it should be appreciated that, since the antenna array 22-2A has 
plural antenna elements, there are a corresponding plurality of channels through which 
wavefronts are received, and accordingly there could also be a separate channel impulse 

10 response or separate channel estimate for each of the plural channels. But by storing 
the complex samples in antenna signal matrix 80 in the manner aforedescribed, and by 
concurrently finding the time of arrival (TO A) and channel coefficients over the entire 
antenna signal matrix 80, the channel estimate generator 60-2A provides a channel 
estimate which encompasses all channels for all antenna elements and for this reason is 

15 known as a "composite" channel estimate. 

[00092] The composite channel estimate comprises, as mentioned before, a time of 
arrival (TO A) and channel coefficient for each arriving wavefront in the sampling 
window (e.g., a channel coefficient mapped to a time of arrival (TOA)). Therefore, the 
channel estimate may comprise a set (of one or more) pairs of data, each pair including 
20 a time of arrival (TOA) and channel coefficient. The task for correlator 50-2A is thus 
to locate a value or "tone" in antenna signal matrix 80 that best corresponds to an 
arriving wavefront, e.g., to locate a value or tone for each arriving wavefront in the 
sampling window. 

[00093] The task of locating a value or "tone" in antenna signal matrix 80 that best 
25 corresponds to an arriving wavefront can be accomplished by various techniques, 

including both parametric and non-parametric techniques. A Fast Fourier Transform 
(FFT) technique as discussed below is just one representative and illustrative example 
non-parametric type of correlator which can be utilized. 

[00094] Fig. 7 depicts example basic steps performed by an example correlator 50-2A 
30 and correlator output analyzer 62-2A in conjunction with the Fast Fourier Transform 
(FFT) calculation. As step 7-1, the correlator 50-2A of Fig. 2A calculates Expression 1. 
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Y(n,t) = FFT(n,X(:,t)) Expression 1 



In Expression 1, t is the sampling window time index; X(:,t) is the complex antenna 
matrix (with the colon representing all antenna indexes for one sampling window 
time index); and, n is the frequency index. Each FFT calculation is thus a one 
5 dimensional FFT calculation on the baseband signal, and corresponds to a specific 
direction of arrival (as depicted by the frequency index) and set of antenna weights 
which, in practice, are the FFT weights. 

[00095] The output of correlator 50-2 A, i.e., the Y(n,t) values computed using 
Expression 1, are stored as correlator output values. The correlator output values can 
10 be stored, for example, in the correlator output value memory 52-2A of Fig. 2A. 

[00096] The correlator output analyzer 62-2A of channel estimate (CE) generator 60- 
2A searches the correlator output values and (as step 7-2) determines therefrom a 
maximum absolute value |Y(n,t)|max- This maximum absolute value |Y(n,t)|max is 
utiUzed by correlator output analyzer 62-2A to determine both the direction of arrival 

15 (DO A) and time of arrival (TOA) for an arriving wavefront seen in the sampling 

window. In particular, as step 7-3 correlator output analyzer 62-2A chooses a sampling 
window time index t_max at which |Y(n,t)|max occurs to be the time of arrival of the 
arriving wavefront. In addition, as step 7-4 correlator output analyzer 62-2A chooses 
the frequency index n_max at which |Y(n,t)|niax occurs to represent the direction of 

20 arrival (DO A) of the arriving wavefront. The frequency index corresponds to a 

direction of arrival (e.g., 0). An amplitude for the arriving wavefront is determined as 
correlator output analyzer 62-2A divides |Y(n,t)|max by the number of antennas 
comprising the antenna array (as step 7-5). 

[00097] Expression 1 and the steps of Fig. 7 represent a generic non-parametric FFT 
25 calculation. In a CDMA-specific situation which utilizes a coding generator (such as 
coding generator 30 of Fig. 1), a comparable FFT calculation can be made using a 
refinement of Expression 1 which appears as Expression 2. 



Y(n,t) = Z Cj *FFT(n,X(:,t)), j = 1,K 



Expression 2 



22 

Expression 2 is understood from Expression 1, it being further mentioned that Cj is a 
coding sequence symbol value j; and K is a length of the coding sequence. 

[00098] As a result of operation of joint searcher and channel estimator 24-2A, an 
accurate channel estimate can be provided to the detector as a spatial signature. The 
5 spatial signature includes the time of arrival (TO A), as well as the direction of arrival 
(DOA) and amplitude. As explained below, the channel coefficient (CC) for each 
wavefront is derived from the direction of arrival (DOA) and amplitude. The time of 
arrival (TO A) and channel coefficient (CC) are applied to the detector as represented by 
lines 66-2 A and 68-2 A, respectively, in Fig. 2A. 

10 [00099] As mentioned above, the channel coefficient (CC) for each wavefront is 
derived from the direction of arrival (DOA) and amplitude. Recall that at step 7-4 
correlator output analyzer 62-2A chose the frequency index n_max at which |Y(n,t)|max 
occurs to represent the direction of arrival (DOA) of the arriving wavefront, with the 
chosen frequency index corresponding to a direction of arrival (e.g., 9). The channel 

15 impulse response vector (i.e., array propagation vector) x is therefore generated by 
detector interface 64-2A in accordance with Expression 3 (for identical isotropic 
antenna elements). 

[000100] X =[l,e^^^*^^"®>, e^'^^*'^^°®>, ... e«^^*<^-i>^^°®>]*C Expression 3 

[000101] In Expression 3, j is the conventional imaginary notation; k = 2*7rX; d is a 
20 spacing distance between elements of the antenna array; X is the wavelength of the 
received/transmitted electro-magnetic signal: (f*X =c) and, K is the antenna element 
index (illustrated as antenna numbers Al, A2, A3, A4 in FIG. 9 A, for example) In 
Expression 3, C is a complex constant in which |C| = |FFT_max [/number of antennas; 
the argument of C, i.e., arg(C) = arg(FFr_max), wherein |FFT_max is the FFT value 
25 computed at step 7-1 of Fig. 7. 

[000102] In the foregoing description it is the role of channel estimate (CE) 
generator 60-2A, and particularly detector interface 64-2A, to generate both a time of 
arrival (TOA) and a channel coefficient (CC), the channel coefficient being derived 
from the direction of arrival, e.g., as above described in conjunction with Expression 3. 
30 In an alternate implementation of this and other embodiments described herein, the 
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detector itself (such as detector 26 illustrated in Fig. 1), upon receiving the time of 
arrival (TO A) and direction of arrival (DO A) for each arriving wavefront, may have the 
intelligence to compute the channel coefficient for each wavefront from the 
corresponding direction of arrival (DOA) information. In such case, the time of arrival 
5 and direction of arrival are output by detector interface 64 to the detector. 

[000103] Thus, considering the aspects above discussed, the joint searcher and 
channel estimator 24-2A looks in a discrete number of possible directions of arrival, 
and picks the direction of arrival with the highest correlation (highest absolute value). 
A comparative operational evaluation was preformed to illustrate the efficacy of a joint 

10 searcher and channel estimator such as joint searcher and channel estimator of 24-2A of 
Fig. 2A. A first scenario of the comparative operational evaluation involved a 
conventional searcher functioned essentially in prior art fashion for a sampling window. 
In so doing, with respect to each antenna for the sampling window the conventional 
search merely picked the time (e.g., chip) which had the greatest absolute value. In 

15 other words, the signals from each antenna were processed separately. A second 

scenario of the comparative operational evaluation was performed in the manner above 
described with respect to the joint searcher and channel estimator 24-2A of Fig. 2 and 
Expression 1 . The same signal was applied in both scenarios to an antenna array 
having eight antenna elements. The length of the sampling window for both scenarios 

20 was twenty chips, and a coding sequence of { 1 } was utilized (e.g., only one of the chips 
contained the signal, the remainder of the chips contained complex white noise). 

[000104] Fig. 8 A illustrates the first scenario which utilized the conventional 
searcher. In contrast. Fig. 8B illustrates the spatial joint searcher and channel estimator 
24-2A of Fig. 2A utihzed for the second scenario. The superiority of the second 

25 scenario (and thus the spatial joint searcher and channel estimator) is evident by a 

comparison of Fig. 8 A and Fig. 8B, due to the higher SNR for the signal of interest in 
Fig. 8B . In the second scenario, it is much easier to pick out the tone or value for the 
arriving wavefront. For the second scenario, Fig. 8C(1) shows the absolute value of the 
complex channel impulse response taps; Fig. 8C(2) shows phase errors of the complex 

30 channel impulse response taps; and Fig. 8C(3) shows the detected time of arrival. 

[000105] Whereas the joint searcher and channel estimator of Fig. 2A includes a non- 
parametric type matrix analyzer, e.g., a correlator (e.g., a filter which performs a Fast 
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Fourier Transform (FFT) calculation), in other example embodiments the matrix 
analyzer of the joint searcher and channel estimator implements parametric techniques. 
As does the Fig. 2A embodiment, the spatial joint searcher and channel estimator 24-2B 
of Fig. 2B (which uses a parametric technique) is shown along with its associated 
5 example antenna array 22-2B. Again by way of example, antenna array 22-2B includes 
four antenna elements 22-2B-1 through 22-2B-4. The signals obtained from the 
antenna elements are each applied to joint searcher and channel estimator 24-2B, as 
well as to a detector (not illustrated in Fig. 2B). 

[000106] Similar to the earlier described embodiment, joint searcher and channel 
10 estimator 24-2B can comprise an antenna signal matrix handling unit 40-2B, which in 
turn comprises antenna signal matrix generator 42-2B and antenna signal matrix 
memory 44-2B, which function much in the manner previously described. For 
example, the complex baseband values stored in antenna signal matrix memory 44-2B 
can also be conceptualized as matrix 80, and as such has a sampling window time 
15 index. The antenna signal matrix 80 has been previously discussed in conjunction with 
Fig. 6, and is now also discussed with reference to Fig. 9A for sake of expounding the 
joint searcher and channel estimator24-2B of Fig. 2B. 

[000107] The joint searcher and channel estimator 24-2B further comprises a matrix 
analyzer, e.g., parametric estimator 51-2B, which utilizes a parametric technique. In 

20 addition, in similar manner as the preceding embodiment, joint searcher and channel 
estimator 24-2B comprises a channel estimate generator 60-2B which has parametric 
estimation output vector analyzer 62-2B and a demodulator interface 64-2B. Basic 
steps performed by parametric estimator 51-2B and parametric estimation output vector 
analyzer 62-2B of the joint searcher and channel estimator 24-2B of Fig. 2B are 

25 illustrated in Fig. 10. 

[000108] For each sampling window time index of the antenna signal matrix 80, as step 
10-1 the parametric estimator 51-2B estimates, e.g. two parameters at each time instant: 
a spatial frequency parameter parameter and a spatial amplitude parameter. The spatial 
frequency parameter estimates the frequency the incident waves creates when arriving 
30 at the ULA. The spatial amplitude parameter estimates the amplitude of this frequency. 
The spatial frequency parameter and spatial amplitude parameter are considered to be a 
parameter pair and in Fig. 9B, they are illustrated as one parameter per sample along 



25 

the sampling time index. The parameters can be calculated by an appropriate strategy 
or goal criteria, e.g., by a minimum mean square error technique (MMSE). 



[000109] As step 10-2, parametric estimation output vector analyzer 62-2B finds certain 
"qualifying" values in parametric estimation output vector, i.e. high or maximum values 
5 of the spatial amplitude parameter. The qualifying values can be, for example, values 
whose absolute values are sufficiently high or are a maximum. Each qualifying value 
of parametric estimation output vector 90 can correspond to an arriving wavefront for 
the sampling window. 

[0001 10] For each qualifying value, as step 10-3 the parametric output estimation 
10 vector analyzer 62-2B chooses a time of arrival (TOA) as corresponding to the 

sampling window time index t for the qualifying value, e.g., the time index at which the 
maximum/qualifying absolute value of the parametric estimation output vector occurs. 

[0001 1 1] Similarly, for each qualifying value, as step 10-4 the analyzer 62-2B chooses 
a direction of arrival (DOA) as the spatial frequency parameter value at the time of 
15 arrival, decided in 10-3. 

[0001 12] As step 10-5, the parametric estimation output vector analyzer 62-2B 
determines the amplitude as the value of the spatial amplitude value divided by the 
number of antenna elements in the array. 

[0001 13] The joint searcher and channel estimator 24-2B thus looks for an 
20 optimum direction, and prepares a channel estimate which can be provided to the 
detector as a spatial signature. The spatial signature includes the direction of arrival 
(DOA) and amplitude. The channel coefficient (CC) for each wavefront is derived 
from the direction of arrival (DOA) and amplitude in the manner explained above with 
reference to Expression 3. The time of arrival (TOA) and channel coefficient (CC) are 
25 applied to the detector as represented by lines 66-2B and 68-2B, respectively, in Fig. 
2B. 

[0001 14] It should be understood from the foregoing that information indicative of 
more than one incident wavefront may be seen in a sampling window. For example, 
with reference to the parametric estimation output vector 90 of Fig. 9B, the parametric 
30 estimation output vector analyzer 62-2B may see other high numbers and for each of 
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those high numbers which qualify, an arriving wavefront may be ascertained. For 
example, if there were two high numbers, then the channel impulse response may 
reflect two arriving wavefronts. For each of the two arriving wavefronts the joint 
searcher and channel estimator would pick out both a time of arrival (TO A) and 
5 direction of arrival (DOA), as well as amplitude, which are mapped to two different 
channel coefficients, with these two different channel coefficients forming part of the 
channel estimate. 



[0001 15] Fig. 4 showed that a wavefront individually reached each of four example 
antenna elements of an antenna array, providing a different antenna output (complex 

10 baseband signal) for each antenna element. For example, the output of antenna 
element 22-1 with the complex vector ai.i (and phase Om); the output of antenna 
element 22-2 is the complex vector ai.2 (and phase 81-2), and so forth. The linear 
combination, of the complex antenna baseband signal and the antenna weight vectors 
Wi have the effect of a summation, or coherent combination in the time and space 

15 domain, shown as summation function 100 in Fig. 12. 

[0001 16] The coherent combination facilitated by the antenna weight vectors Wj is 
illustrated in Fig. 11. In the example case of the four antenna elements shown in Fig. 
12, the effect of weight belonging to antenna index 2, here denoted asW2 is to rotate the 
output of antenna element 22-2 so that its phase 61.2 lines up with the phase 0 1-1 of the 

20 output of antenna element 22-1, in the manner shown in Fig. 1 1 . Similarly, the effect of 
weight W3 is to rotate the output of antenna element 22-3 so that its phase 0i_3 lines up 
with the phase Oi.j of the output of antenna element 22-1. The effect of weight W4 is 
to rotate the output of antenna element 22-4 so that its phase 81.4 lines up with the phase 
0i_i of the output of antenna element 22-1. For simplicity, Fig. 1 1 ignores noise 

25 considerations, which tend to make the resultant vector less than straight. Note that in 
the preciding paragraph, the weight vectors are denoted with Wi, where i denotes the 
antenna index of the weight vector W, which is denoted without index. 

[0001 17] In the spatial joint searcher and channel estimators, the SINR for finding 
channel taps (peaks) should be proportional to the number of antenna elements 
30 comprising the array. The operation of the spatial joint searcher and channel estimators 
can be adapted to take into consideration channel variations over time, e.g., spatial 
variations in the environment (e.g., in the sending and receiving antennas). 
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[0001 18] The non-parametric FFT-type correlator and the parametric techniques 
illustrated above, e.g., by Fig. 2A and Fig. 2B, respectively, are only two example 
techniques for finding the values or "tones" in antenna signal matrix 80 which are 
associated with arriving wavefronts. Other parametric approaches are described by or 
5 understood from Stocia, Petre and Moses, Randolph, Introduction To Spectral Analysis, 
ISBN-013-258419-0, Prentice Hall, which is incorporated by reference in its entirety, 
particularly Chapter 4 thereof. 

[0001 19] The spatial joint searcher and channel estimator and techniques of 
operation thereof as described above are suitable for any receiver unit which has plural 
10 receiving antennas. Thus, the spatial joint searcher and channel estimator is particularly 
well suited for, but not limited to, a base station which has plural antennas. Also 
encompasses are mobile terminals which have plural antennas. 

[000120] TEMPORAL JOINT SEARCHERS/ESTIMATORS 

[000121] In other embodiments, the joint searcher and channel estimator 
15 contemporaneously processes the signals received at an antenna element from plural, 
successive sets of pilot data (each set of pilot data being received in its own sampling 
window) in order to determine both time of arrival and the channel coefficient. In so 
doing, the joint searcher and channel estimator takes into consideration a doppler shift 
or frequency shift (the terms "doppler shift" and "frequency shift" being used 
20 interchangeably in conjunction with the description of the temporal joint searcher and 
channel estimator). The frequency shift is primarily attributable to a doppler shift, but 
can also include a frequency shift in the transmitter and receiver oscillators. For 
simplification, such frequency shifts are hereinafter referred to as "doppler shifts" or 
"doppler frequency shifts " . 

25 [000122] The doppler shift can be occasioned by movement such as relative 

movement of one of the transmitter and the receiver (for example, by movement of a 
mobile terminal), or movement of a signal path-affecting object or structure in the 
surroundings (which can cause a doppler shift even for a fixed transmitter and fixed 
receiver). 



[000123] In providing the channel estimate, the joint searcher and channel 
estimator essentially concurrently considers plural signals (e.g., plural sets of pilot data) 
received by the antenna element. The joint searcher and channel estimator applies the 
channel coefficient and the time of arrival to a detector which provides, e.g., a symbol 
5 estimate. 

[000124] In these embodiments, the joint searcher and channel estimator is 
essentially a two dimensional unit, with a second dimension being a temporal 
dimension imparted by the time intervals at which the successive sets of pilot data 
arrive. This temporal dimension, which involves essentially simultaneous and 
10 concurrent processing together of signals received at the antenna element from each of 
the plural sets of pilot data, bestows on these embodiments of the joint searcher and 
channel estimator the distinction of being a "temporal" joint searcher and channel 
estimator. 

[000125] The temporal joint searcher and channel estimator may take differing 
15 embodiments and have differing implementations. In one example, illustrative 
embodiment the temporal joint searcher and channel estimator includes a non- 
parametric type correlator (e.g., a correlator which performs a Fast Fourier Transform 
(FFT) calculation). In another example, illustrative embodiment the temporal joint 
searcher and channel estimator utilizes a parametric approach. 

20 [000126] Fig. 13A illustrates one example embodiment of a spatial joint searcher and 
channel estimator 24-13A which uses a non-parametric technique for determining time 
of arrival and channel estimate, as well as an associated example antenna array 22- 13 A. 
In the example of Fig. 13 A, the antenna array 22-1 3A is shown as having one antenna 
element 22-13A-1. As explained hereinafter, complex baseband signals obtained from 

25 the same antenna element (e.g., antenna element 22-13A-1) upon receipt of each of 

successive sets of pilot data (as hereinafter described) are each applied to joint searcher 
and channel estimator 24-13A, as well as to a detector (not illustrated in Fig. 13 A). 

[000127] The joint searcher and channel estimator 24-13A comprises an antenna signal 
matrix handling unit 40- 13 A. In one particular example manifestation, antenna signal 
30 matrix handling unit 40-13A includes antenna signal matrix generator 42-1 3A and 
antenna signal matrix memory 44- 13 A. A matrix analyzer, which for the non- 
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parametric technique of Fig. 2A can be correlator 50- 13 A, operates on complex values 
stored in antenna signal matrix memory 44- 13 A. The correlator 50-13A preferably 
comprises a filter. The correlator 50-13A generates certain output values, which may 
be stored, e.g., in correlator output value memory 52- 13 A. The joint searcher and 
5 channel estimator 24-1 3A further comprises a channel estimate (CE) generator 60- 13 A. 
In the illustrated example embodiment, the chanouel estimate (CE) generator 60-1 3 A 
comprises a correlator output analyzer 62-13A and a detector interface 64- 13 A. The 
detector interface 64-13A generates, for each wavefront, a channel estimate which 
includes both a time of arrival (TO A) and a channel coefficient (CC). In Fig. 13 A, the 
10 time of arrival and channel coefficient output by detector interface 64-1 3A are applied 
to the detector on lines 66-13A and 68- 13 A, respectively. 

[000128] As shown in Fig. 14, the temporal joint searcher and channel estimators such 
as joint searcher and channel estimator 24-13A of Fig. 13A watch the channel response 
from an antenna (e.g., antenna 22-13-1) for sets of pilot data which are interpersed or 

15 otherwise transmitted with other data (e.g., user data). For sake of simplicity, it is 
assumed that each set of pilot data is received in a separate sampling window. Such 
need not be the case, however, as differing sets of pilot data can be received 
simultaneously if the different streams are, e.g., code multiplexed. Merely as an 
illustrative example, Fig. 14 shows four sets of pilot data, i.e., pilot sets Tl - T4, 

20 interspersed with user data and received at unique global times (as indicated by the "T" 
axis in Fig. 4). 

[000129] Each set of pilot data is typically in a different frame from another set of pilot 
data. For example, pilot set Tl may be in frame 1 ; pilot set T2 may be in frame 1 1 ; 
pilot set T3 may be in frame 21 ; etc. "Frame transmission interval" refers to the time 
25 between two successive frames which contain pilot data. The time between two 

successive frames which contain pilot data is typically specified by a standard or other 
specification. 

[000130] Fig. 14 thus reflects the typical periodic transmission of pilot data by the 
transmitter source, and also the expected receipt of repetitions of the pilot data at the 
30 receiver at successive intervals. In view of factors such as relative motion of the 

transmitter and receiver, the successive intervals between differing sets of pilot data are 
not necessarily constant. 
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[000131] As further shown in Fig. 14, an antenna matrix handling unit (such as antenna 
matrix handling unit 40-13A of the Fig. 13A embodiment) samples the signals received 
by the antenna element for each of the successive sets of pilot data, i.e., for pilot sets Tl 
- T4. Using the sampled signals, antenna signal matrix generator 42-13A generates an 
5 antenna signal matrix such as antenna signal matrix 110 illustrated in Fig. 14. The 
antenna signal matrix 110 may be stored in any convenient fashion, such as antenna 
matrix memory 44- 13 A. 

[000132] The antenna signal matrix 1 10 is a two dimensional functionally dependent 
matrix. In other words, complex samples are stored in antenna signal matrix 1 10 as a 

10 function of two different indexes. For the antenna signal matrix 110 shown in Fig. 14, 
a first index is a sampling window time index, illustrated along the X axis of Fig. 14. 
For embodiments which utilize spreading codes or similar codes, the first index may be, 
for example, a chip index. Thus, the sampling window time index points to a time in 
the sampling window relative to a start of the respective sampling window. In the 

15 antenna signal matrix 1 10 of Fig. 14, a second index, shown along the Y axis, is a pilot 
set index (which serves as a dimensional differentiation index). The pilot set index 
indicates which one of the sets of pilot data the sample was obtained. In other words, a 
pilot set index = Tl indicates that the sample was obtained from pilot set Tl; a pilot set 
index = T2 indicates that the sample was obtained from pilot set T2; and so on as 

20 depicted by the arrows which connect the matrix 110 with the received signal depiction 
with its illustrative successive sets of pilot data. As can be seen, the pilot set index 
points to a different row of the antenna signal matrix 110, each row being associated 
with a different set of pilot data. 

[000133] Fig. 14 shows four rows in antenna signal matrix 110 for consistency with the 
25 illustrated example wherein the antenna signal matrix encompasses four successive sets 
of pilot data. The number of sets of pilot data subsumed in a given antenna signal 
matrix, and thus the maximum value of the pilot set index, can vary from receiver to 
receiver, so that the present example's choice of four sets of pilot data is only 
illustrative for sake of example. In general, the choice of the number of sets of pilot 
30 data to be apprehended simultaneously by a temporal joint searcher and channel 

estimator depends on how quickly the doppler is expected to change. The number of 
taps/incident waves depends on the multipath. In other words in an open space we have 
one direct path and thus only one channel/tap coefficient in the channel impulse 
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response. 

[000134] The antenna signal matrix 1 10 is also conceptualized as storing 
"dimensionally differentiated" signals acquired from a single antenna element of the 
antenna array. For the temporal joint searcher and channel estimator, wherein the 
5 antenna structure comprises an antenna which provides signals for each of successive 
sets of pilot data received at separated time intervals, the signals acquired by the 
antenna are dimensionally differentiated with regard to a temporal or time dimension. 
For example, the signals acquired by the antenna are dimensionally differentiated by 
being acquired in differing frame transmission intervals. 

10 [000135] For sake of simplicity, the complex values stored in antenna signal matrix 

110, including the complex values obtained from the antennas, are not illustrated in Fig. 
14. Such complex values would be illustrated in a third dimension, e.g., out of the 
plane of Fig. 14. The antenna signal matrix 110 includes both complex white noise and 
(for the sake of the present illustration) a complex sample for at least one wavefront 

15 (planar or other known shape). The wavefronts have known phase (temporal, non- 
coherent detection), and are modulated code sequences. 

[000136] The complex values stored for each column of the antenna signal matrix 
1 10 of Fig. 14 can be conceptualized as a dimensional receptivity vector. That is, a 
dimensional receptivity vector is formed with the complex values taken with respect to 

20 a same single sampling window time index for each of the sets of pilot signals included 
in the sampling window (e.g., for sets Tl - T4 in Fig. 14). Each element taken from a 
unique row of antenna signal matrix 110 has a different phase in the manner of the 
differing 0 values illustrated in Fig. 5. As received by the differing antenna elements, 
for the temporal joint searcher and channel estimator the change in phase over time is 

25 the doppler frequency for the dimensional receptivity vector. The phase rotation speed, 
or frequency, of the dimensional receptivity vector, for the sampling window time 
instance can be interpreted as the doppler shift (DS). Thus, each dimensional 
receptivity vector corresponds to a separate doppler shift frequency. There are plural 
possible frequencies for the dimensional receptivity vector, each of the plural possible 

30 frequencies corresponding to a different possible doppler shift for a wavefront. For the 
non-parametric techniques herein employed, the plural possible frequencies can be a 



32 



continuous range of frequencies. For sake of differentiating the plural possible 
frequencies, the plural possible frequencies are each represented by a frequency index. 

[000137] For the temporal joint searcher and channel estimator, the channel 
estimate comprises, as mentioned before, a time of arrival (TO A) and doppler shift for 
5 each arriving wavefront in the sampling window (e.g., a channel coefficient mapped to 
a doppler shift). Therefore, the channel estimate may comprise a set (of one or more) 
pairs of data, each pair including a time of arrival (TO A) and a channel coefficient. 
The task for the temporal joint searcher and channel estimator is thus to locate a value 
or "tone" in antenna signal matrix 110 that best corresponds to an arriving wavefront, 

10 e.g., to locate a value or tone for each arriving wavefront in the sampling window. This 
task of locating a value or "tone" in antenna signal matrix 110 that best corresponds to 
an arriving wavefront can be accomplished by various techniques, including both 
parametric and non-parametric techniques. A Fast Fourier Transform (FFT) technique 
as discussed below is just one representative and illustrative example non-parametric 

15 type of correlator which can be utilized. 

[000138] Fig. 15 depicts example basic steps performed by an example correlator 
50-1 3A and correlator output analyzer 62-13A in conjunction with the Fast Fourier 
Transform (FFT) calculation. As step 15-1, the correlator 50-13A of Fig. 13A 
calculates Expression 5. 

20 Y(n,t) = FFT(n,X(n,t)) Expression 5 

wherein t is the sampling window time index; X(n,t) is the complex antenna matrix; 
and, n is the doppler frequency index. Each FFT calculation is thus a one dimensional 
FFT calculation on the baseband signal, and corresponds to a specific doppler shift 
frequency. 

25 [000139] The output of correlator 50-13A, i.e., the Y(n,t) values computed using 

Expression 1, are stored as correlator output values. The correlator output values can 
be stored, for example, in the correlator output value memory 52-13A of Fig. 13 A. 

[000140] The correlator output analyzer 62-1 3A of channel estimate (CE) generator 
60-13A searches the correlator output values and (as step 15-2) determines therefrom a 
30 maximum absolute value |Y(n,t)|max- This maximum absolute value |Y(n,t)|niax is 
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utilized by correlator output analyzer 62-1 3A to determine both the doppler shift (DS) 
and time of arrival (TO A) for an arriving wavefront. In particular, as step 15-3 the 
correlator output analyzer 62-1 3A chooses a sampling window time index t_max at 
which |Y(n,t)|max occurs to be the time of arrival of the arriving wavefront. In addition, 
5 as step 15-4 the correlator output analyzer 62-13A chooses the doppler index n_max at 
which |Y(n,t)|niax occurs to determine the doppler shift (DS) of the arriving wavefront. 
An amplitude for the arriving wavefront is determined as correlator output analyzer 62- 
13A divides |Y(n,t)|niax by the number of sets of pilot data comprising the antenna signal 
matrix (as step 15-5). 

10 [000141] Expression 5 and the steps of Fig. 15 represent a generic FFT calculation. 
In a CDMA-specific situation which utilizes a coding generator (such as coding 
generator 30 of Fig. 1), a comparable FFT calculation can be made using a refinement 
of Expression 5 such as that which appears as Expression 2, previously discussed, but 
applied for the temporal joint searcher and channel estimator rather than for the spatial 

15 joint searcher and channel estimator. 

[000142] As a result of operation of joint searcher and channel estimator 24-13A, 
an accurate channel estimate can be provided to the detector as a temporal signature. 
For each wavefront, the temporal signature includes the time of arrival (TO A) mapped 
to a doppler (frequency) shift. As explained below, the channel coefficient (CC) for 
20 each time of arrival and wavefront is derived from the doppler frequency shift. The 
time of arrival (TO A) and channel coefficient (CC) are applied to the detector as 
represented by Unes 66-13A and 68-13A, respectively, in Fig. 13A. 

[000143] As mentioned above, the channel coefficient (CC) for each wavefront is 
derived from the doppler frequency shift (DS). Recall that at step 15-4 the correlator 
25 output analyzer 62-2B chose the frequency index n_max at which |Y(n,t)|niax occurs to 
represent the doppler shift frequency (DSF) of the arriving wavefront, with the chosen 
frequency index corresponding to a doppler shift (e.g., 9', i.e., the derivative of 9). The 
channel impulse response vector (i.e., array propagation vector) x is therefore generated 
by detector interface 64-2B in accordance with Expression 6. 

30 [000144] C[ei'^ e*'"^' "\ d'^^ ^ . . .e*'"^ ^ "^] Expression 6 
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[000145] In Expression 6, C is the amplitude of the wavefront, f is the frequency of 
the signal (including doppler shift); T is the period time between two pilot 
symbols/sequences (which are assumed to be periodical, in analogy to the uniform 
array of the spatial embodiment), and H is a complex value of the signal at the first pilot 
5 symbol/sequence, H being arg(FFT max). For sake of simpUcity, noise has been 
excluded from Expression 6, and C is assumed to be constant within the time TN. 

[000146] In the foregoing description it is the role of channel estimate (CE) 
generator 60-2A, and particularly detector interface 64-2A, to generate both a time of 
arrival (TO A) and a channel coefficient (CC), the channel coefficient being derived 

10 from the doppler shift, e.g., as above described in conjunction with Expression 6. In an 
altemate implementation of this and other embodiments described herein, the detector 
itself (such as detector 26 illustrated in Fig. 1), upon receiving the time of arrival 
(TO A) and doppler shift (DS) for each arriving wavefront, may have the intelligence to 
compute the channel coefficient for each wavefront from the corresponding direction of 

15 arrival (DOA) information. In such case, the time of arrival and direction of arrival are 
output by detector interface 64-13A to the detector. 

[000147] Thus, the joint searcher and channel estimator 24-13A looks at a discrete 
number of possible doppler frequency shifts, and picks the doppler frequency with the 
highest correlation (highest absolute value). 

20 [000148] Whereas the joint searcher and channel estimator of Fig. 13A includes a non- 
parametric correlator (e.g., a filter) which performs a Fast Fourier Transform (FFT) 
calculation, in other example embodiments the temporal joint searcher and channel 
estimator implements parametric techniques. As does the Fig. 13A embodiment, the 
spatial joint searcher and channel estimator 24-13B of Fig. 13B is shown along with its 

25 associated example antenna array 22-13B comprising an antenna element 22-13B-1 
which receives the successive sets of pilot data in the manner of Fig. 14. 

[000149] Similarly to the earlier described embodiment, joint searcher and channel 
estimator 24-13B can comprise an antenna signal matrix handling unit 40-13B, which 
in turn comprises antenna signal matrix generator 42-13B and antenna signal matrix 
30 memory 44-13B, which function much in the manner previously described. For 

example, the complex baseband values stored in antenna signal matrix memory 44-13B 



can also be conceptualized as matrix 110, and as such has a sampling window time 
index. The antenna signal matrix 110 has been previously discussed in conjunction 
with Fig. 14, and is now also discussed with reference to Fig. 16A for sake of 
expounding the joint searcher and channel estimator 24-13B of Fig. 13B. 

5 [000150] The joint searcher and channel estimator 24-13B further comprises the 

parametric estimator 51-13B which outputs a parametric output estimation vector for 
storage in memory 52-13B. In addition, in similar manner as the preceding 
embodiment, joint searcher and channel estimator 24-13B comprises a channel estimate 
generator 60-13B which has parametric output estimation vector analyzer 62-13B and a 
10 demodulator interface 64- 1 3B . Basic steps performed by parametric estimator 5 1 - 1 3B 
and parametric output estimation vector analyzer 62-13B of the joint searcher and 
channel estimator 24-13B of Fig. 13B are illustrated in Fig. 17. 

[000151] For each sampling window time index of the antenna signal matrix 1 10. As 
step 17-1, the parametric estimator 51-13B estimates, e.g., two parameters at each time 

15 instant: a temporal frequency parameter parameter and a temporal amplitude 

parameter. The temporal frequency parameter estimates the frequency the incident 
waves creates when arriving at the antenna for the consecutive pilot symbols. The 
temporal amplitude parameter estimates the amplitude of this frequency. The temporal 
frequency parameter and temporal amplitude parameter are considered to be a 

20 parameter pair and in Fig. 16B , they are illustrated as one parameter per sample along 
the sampling time index. 

[000152] As step 17-2 performed by joint searcher and channel estimator 24-13B, 
analyzer 62-13B finds certain "qualifying" values in parametric output estimation 
vector 120, i.e. maximum value of temporal amplitude vector. Each qualifying value of 
25 parametric output estimation vector 120 can correspond to an arriving wavefront for the 
sampling window. 

[000153] For each qualifying value, as step 17-3 the parametric output estimation 
vector analyzer 62-13B chooses a time of arrival (TO A) as corresponding to the 
sampling window time index t for the qualifying value,, e.g., the time index at which the 
30 maximum/qualifying absolute value of the parametric estimation output vector occurs 
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[000154] Similarly, for each qualifying value, as step 17-4 the parametric output 
estimation vector analyzer 62-13B chooses a a doppler shift frequency (DS) as the 
temporal frequency parameter value at the time of arrival decided in 17-3. 

[000155] As step 17-5, the parametric estimation output vector analyzer 62-1 3B 
5 determines the amplitude as being the maximum/qualifying absolute value divided by 
the number of pilot data sets in the series. 

[000156] The joint searcher and channel estimator 24-13B thus looks for an 
optimum doppler (shift) frequency, and prepares a channel estimate which can be 
provided to the detector as a temporal signature. The temporal signature includes the 
10 time of arrival (TO A), as well as the doppler shift frequency (DSF) and amplitude. The 
channel coefficient (CC) for each time of arrival and wavefront is derived from the 
doppler shift (DS) in the manner described above with reference to Expression 6. The 
time of arrival (TO A) and channel coefficient (CC) are applied to the detector as 
represented by lines 66-13B and 68-13B, respectively, in Fig. 13B. 

15 [000157] It should be understood from the foregoing that information indicative of 
more than one incident wavefront may be seen in a sampling window. For example, 
with reference to the parametric output estimation vector 120 of Fig. 16B, the 
parametric output estimation vector analyzer 62-13B may see other high numbers and 
for each of those high numbers which qualify, an arriving wavefront may be 

20 ascertained. For example, if there were two high numbers, then the channel impulse 
response may reflect two arriving wavefronts. For each of the two arriving wavefronts 
the joint searcher and channel estimator would pick out both a time of arrival (TO A) 
and doppler shift frequency (DSF), as well as amplitude, which are mapped to two 
different channel coefficients, with these two different channel coefficients forming 

25 part of the channel estimate for the channel impulse response. 

[000158] The operation of the temporal searcher and channel estimator has been 
described above for one antenna element of an antenna array 22. It should be 
understood that the antenna array 22 may comprise plural antenna elements, and that 
the operations described above may be performed separately with respect to one or 
30 more antenna elements of the array. Moreover, as described later, principles of the 
foregoing operation may be performed in a combined manner with respect to plural 
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antennas of the antenna array. 

[000159] The temporal joint searcher and channel estimator and techniques of 
operation thereof as described above is particularly well suited for, but not limited to, a 
receiver unit which has only one antenna element, e.g., a mobile terminal with only one 
5 antenna. As indicated above, however, the temporal joint searcher and channel 

estimation techniques can be utilized separately but in parallel by plural antennas for a 
receiver. 

[000160] Consider, for example, the situation reflected in Fig. 1 1 in which the output of 
antenna element 22-13A-1 (or 22-13B-1) for a pilot data set Tl with the complex vector 

10 ai.i (and phase Gm); the output of the same antenna element for a pilot data set T2 with 
the complex vector ai.2 (and phase 9i.2), and so forth. In this situation, the linear 
combinationof the complex antenna baseband signal and the doppler weight vectors Wi 
also has the effect of a summation, or coherent combination in the time domain, shown 
as summation function 100 in Fig. 12. By adding these complex vectors coherently, the 

15 temporal joint searcher and channel estimator increases the performance of the search 
and the channel estimate. 

[000161] In situations in which there is no doppler shift (e.g., the mobile terminal 
stands still or moves in a radial direction relative to the base station), the doppler shift 
frequency may be zero. In such cases the pilot data of the arriving wavefront(s) have 

20 essentially the same complex values. The situation of no doppler shift is just one 
special case of the generic operation of the temporal joint searcher and channel 
estimator described above. When the mobile starts to move, a doppler shift may occur, 
the temporal joint searcher and channel estimators obtains the doppler shift frequency, 
and thereby enhances the channel estimate. The channel estimate is enhanced by 

25 considering the doppler shift, regardless of the magnitude of the doppler shift. 

[000162] The non-parametric FFT-type correlator and the parametric estimator 
techniques illustrated above, e.g., by Fig. 13A and Fig. 13B, respectively, are only two 
example techniques for finding the values or "tones" in antenna signal matrix 110. 
Other parametric approaches are described by or understood from Stocia, Petre and 
30 Moses, Randolph, Introduction To Spectral Analysis, ISBN-013-258419-0, Prentice 
Hall, which is incorporated by reference in its entirety, particularly Chapter 4 thereof. 
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[000163] SPATIAL-TEMPORAL JOINT SEARCHERS/ESTIMATORS 



[000164] In some further embodiments, which combine features from both the 
spatial and temporal embodiments discussed above, plural antenna elements of an 
antenna array provide respective plural series of signals for successive sets of pilot data. 
5 The joint searcher and channel estimators of these further embodiments essentially 
concurrently consider the plural series of signals provided by the plural antennas for 
determining both a time of arrival and channel coefficient. 

[000165] By concurrently considering the signals provided by plural antennas, the 
channel estimate takes into consideration direction of arrival in determining the time of 

10 arrival and channel coefficient. By concurrently considering the series of signals 

provided by each antenna, in which each series comprises successive sets of pilot data, 
the channel estimate further takes into consideration a frequency shift which may be a 
doppler shift (occasioned by relative movement of a transmitter and the receiver or of 
an object in a field between the transmitter and receiver). The channel estimate is 

15 performed by considering spatial and temporal domain jointly and concurrently. 

[000166] Since it processes the series of signals from plural antennas, with each 
series comprising successive sets of pilot data, the joint searcher and channel estimator 
is considered a three dimensional unit. A first dimension is with reference to a time 
index of a sampling window, i.e., a sampling window time index. A second dimension 

20 is a spatial dimension imparted by the spacing of the plural antennas of the array. This 
spatial dimension, which involves essentially simultaneous and concurrent processing 
together of signals from the plural antennas for the array in order to determine the time 
of arrival and channel coefficient, bestows on the joint searcher and channel estimator 
the distinction of being a "spatial" joint searcher and channel estimator. A third 

25 dimension is a temporal dimension imparted by the time interval reflected by the 
successive sets of pilot data. This temporal dimension, which involves essentially 
simultaneous and concurrent processing together of signals for each of the successive 
sets of pilot data in order to determine the time of arrival and channel coefficient, 
bestows on the joint searcher and channel estimator the distinction of being a 

30 "temporal" joint searcher and channel estimator. In view of being both a spatial and 
temporal joint searcher and channel estimator, the joint searcher and channel estimator 
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is also referred to as a "combined" spatial/temporal joint searcher and channel 
estimator, or spatio/temporal joint searcher and channel estimator. 



[000167] Concurrent consideration of the plural series of signals may be either in a 
three dimensional essentially concurrent mode or a sequenced mode. The three 
5 dimensional essentially concurrent mode involves a single step determination of the 
time of arrival and channel coefficient by simultaneously considering signals from all 
antennas of the array for all of the plural series. The sequenced mode involves a two 
step determination of the time of arrival and channel coefficient. In the sequenced 
mode, a first step comprises determining a time of arrival and direction of arrival by 

10 concurrently considering the plural signals provided by the plural antennas for a first of 
the plural series. A second step of the sequenced mode comprises refining the 
estimation of the channel coefficient based on doppler shift by concurrently considering 
the elements of the plural series having the direction of arrival determined in the first 
step. This procedure could also be performed the other way around: first determining 

15 the time of arrival and Doppler shift and then refining the channel estimate by 

concurrently considering the elements of the plural series having the Doppler shift 
determined in the first step. 

[000168] Fig. 18A illustrates one example embodiment of a spatio-temporal joint 
searcher and channel estimator 24- 13 A, as well as an associated example antenna array 

20 22- ISA. The antenna array 22- ISA includes, by way of non-limiting example, four 
antenna elements 22-18A-1 through 22-1SA-4. While the antenna elements 22-lSA-l 
through 22-1 8 A-4 are shown as fonning a uniform linear array (ULA), it should be 
understood that antenna configurations other than a uniform linear are possible, and that 
the number of antenna elements in the antenna array may vary (e.g., the number of 

25 antenna elements is not limited to four). After suitable radio frequency processing, the 
signals obtained from the antenna elements are each applied as baseband signals to joint 
searcher and channel estimator 24- ISA, as well as to a detector (not illustrated in Fig. 
ISA). 

[000169] The joint searcher and channel estimator 24- ISA comprises an antenna signal 
30 matrix handling unit 40- ISA. In one particular example manifestation, antenna signal 
matrix handling unit 40- ISA includes antenna signal matrix generator 42- ISA and 
antenna signal matrix memory 44- ISA. A matrix analyzer, which for the non- 
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parametric technique of Fig. 18A C£in be a correlator 50- 18 A, operates on complex 
values stored in antenna signal matrix memory 44- ISA. The correlator 50- ISA 
preferably comprises a filter. The correlator 50- ISA generates certain output values, 
which may be stored, e.g., in correlator output value memory 52- ISA. The joint 
5 searcher and channel estimator 24-1 8A further comprises a channel estimate (CE) 
generator 60- ISA. In the illustrated example embodiment, the channel estimate (CE) 
generator 60- ISA comprises a correlator output analyzer 62- ISA and a detector 
interface 64- ISA. The detector interface 64-1 SA generates, for each wavefront, a 
channel estimate which includes both a time of arrival (TO A) and a channel coefficient 
10 (CC). In Fig. ISA, the time of arrival and channel coefficient output by detector 
interface 64 are applied to the detector on lines 66-1 SA and 6S-1SA, respectively. 

[000170] In the joint searcher and channel estimator 24- ISA of Fig. ISA, for each series 
of sets of pilot data (represented by pilot data sets Tl - T4), the antenna matrix handling 
unit 40- ISA samples the signals from each antenna element. Using the sampled 
15 signals, antenna signal matrix generator 42- ISA generates an antenna signal matrix 
such as antenna signal matrix 130 illustrated in Fig. 19. The antenna signal matrix 130 
may be stored in any convenient fashion, such as antenna matrix memory 44- ISA. 

[000171] The antenna signal matrix 130 is a three dimensional functionally dependent 
matrix. In other words, complex samples are stored in antenna signal matrix 130 as a 
20 function of three different indexes. For the antenna signal matrix 130 shown in Fig. 19, 
a first index is a sampling window time index, illustrated along the X axis of Fig. 19. 
For embodiments which utilize spreading codes or similar codes, the first index may be, 
for example, a chip index. Thus, the sampling window time index points to a time in 
the sampling window relative to a start of the sampling window. 

25 [000172] In the antenna signal matrix 130 of Fig. 19, a second index, shown along the 
Y axis, is an antenna index. The antenna index points to a different row of the antenna 
signal matrix 130, each row being associated with a different antenna element in 
antenna array 22. Fig. 19 shows four rows in antenna signal matrix 130 for consistency 
with the previous examples of an antenna array comprising four antenna elements. It is 

30 reiterated, however, that the number of antennas in an antenna array, and thus the 
number of rows in antenna signal matrix 130 and the maximum value of the antenna 
index, can vary from receiver to receiver, and that the choice of four antenna is only 
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illustrative for sake of example. 

[000173] In the antenna signal matrix 130 of Fig. 19, a third index, shown along the Z 
axis, is a pilot set index. The pilot set index indicates which one of the sets of pilot data 
the sample was obtained. In other words, a pilot set index = Tl indicates that the 
5 sample was obtained from pilot set Tl ; a pilot set index = T2 indicates that the sample 
was obtained from pilot set T2; and so on as depicted by the arrows which connect the 
matrix 110 with the received signal depiction with its illustrative successive sets of pilot 
data. As can be seen, the pilot set index points to a different plane of the antenna signal 
matrix 110, each plane being associated with a different set of pilot data. 

10 [000174] Fig. 19 shows four planes in antenna signal matrix 130 for consistency with 
the illustrated example wherein the antenna signal matrix encompasses four successive 
sets of pilot data. The number of sets of pilot data subsumed in a given antenna signal 
matrix, and thus the maximum value of the pilot set index, can vary from receiver to 
receiver, so that the present example's choice of four sets of pilot data is only 

15 illustrative for sake of example. In general, the choice of the number of sets of pilot 
data to be apprehended simultaneously by a spatio-temporal spatial/temporal joint 
searcher and channel estimator depends on how quickly the doppler is expected to 
change. The number of taps/incident waves depends on the multipath. In other words 
in an open space we have one direct path and thus only one channel/tap coefficient in 

20 the channel impulse response. 

[000175] For sake of simplicity, the complex values stored in antenna signal matrix 
130, including the complex values obtained from the antennas, are not illustrated in Fig. 
19. Such complex values would be illustrated in a fourth dimension. 

[000176] In conjunction with the antenna signal matrix 130 of Fig. 19, and particularly 
25 a WCDMA case in which spacing of the antenna elements in the antenna array is not 
too far apart, the plane wavefront arriving at the antenna array can be considered to 
arrive in the same sampling window time index (or chip index). 

[000177] Assuming that the wavefront arrives at the antenna elements at different 
times (the time differences are small in comparison to the sampling time interval), the 
30 complex values stored for each column of the antenna signal matrix 130 of Fig. 19 have 
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differing phase (e.g., 9) values in each row of the column. For a uniformly spaced array 
of antenna elements, the phase difference is essentially the same between adjacent rows 
of the same column (although noise may be a factor). But whatever the spacing, the 
rate of change of the phase with respect to time (time of travel of the approaching 
5 wavefront) is the phase rotation speed, or frequency, for the vector formed by the 
colunm, as previously explained. This per column frequency can be interpreted as a 
direction of arrival (DOA). There are plural possible frequencies for the columns of 
antenna signal matrix 130, with each of the plural possible frequencies corresponding to 
a possible direction of arrival (DOA) of a wavefront. The plural possible direction of 
10 arrival frequencies are represented by a frequency index "ni". 

[000178] In a similar manner, for each slice of antenna signal matrix 130 along the 
Z" direction the complex values have differing phase (e.g., 9) values. The Z-aligned 
elements of differing "Z" planes of antenna signal matrix 130 have differing phase 
values in view of a possible doppler shift as detected by the different sets of pilot data 

15 as gathered over plural sets of pilot data in a series. The rate of change over time of the 
phase along the Z direction between successive sets of pilot data is a frequency which is 
associated with the doppler shift. There are plural possible frequencies for the Z slices 
of antenna signal matrix 130, with each of the plural possible frequencies 
corresponding to a possible doppler shift (DS) for wavefront. The plural possible 

20 doppler shift frequencies are represented by a frequency index "n2". 

[000179] The channel estimate generator 60-18A (see Fig. ISA) seeks to develop a 
"composite" channel estimate based on the complex values stored in antenna signal 
matrix 130. As mentioned before, since an antenna array such as antenna array 22- ISA 
has plural antenna elements, there are a corresponding plurality of channels through 

25 which wavefronts are received, and accordingly there could also be a separate channel 
impulse response or separate channel estimate for each of the plural channels. But by 
storing the complex samples in antenna signal matrix 130 in the manner aforedescribed, 
and by concurrently finding the time of arrival (TO A) and channel coefficients over the 
entire antenna signal matrix 130, the channel estimate generator 60-1 8A provides a 

30 channel estimate which encompasses all channels for all antenna elements and for this 
reason is known as a "composite" channel estimate. 
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[000180] The composite channel estimate comprises, as mentioned before, a time 
of arrival (TO A) and channel coefficient for each arriving wavefront in the sampling 
window (e.g., a channel coefficient mapped to a time of arrival (TOA)). Therefore, the 
channel estimate may comprise a set (of one or more) pairs of data, each pair including 
5 a time of arrival (TOA) and channel coefficient. The task for correlator 50-18A is thus 
to locate a value or "tone" in antenna signal matrix 130 that best corresponds to an 
arriving wavefront, e.g., to locate a value or tone for each arriving wavefront in the 
sampling window. 

[000181] The task of locating, in an antenna signal matrix such as antenna signal 
10 matrix 130, a value or "tone" that best corresponds to an arriving wavefront can be 
accomplished by various techniques, including both parametric and non-parametric 
techniques. A Fast Fourier Transform (FFT)) performed in a three dimensional 
essentially concurrent mode is discussed below in conjunction with as just one 
representative and illustrative example of a non-parametric technique wherein 
15 correlator 50- 1 8 A is utilized. 

[000182] Fig. 20 depicts example basic steps performed by an example correlator 
50-18A and analyzer 62-18A in conjunction with the Fast Fourier Transform (FFT) 
calculation. In conjunction with Fig. 20, Fig. 21 shows an antenna signal matrix; a 
doppler weight and antenna weight vector; and a non-parametric estimation output 
20 vector for an example embodiment of a spatio-temporal joint searcher and channel 
estimator which operates in a three dimensional essentially concurrent mode. As step 
20-1, the correlator 50-18A of Fig. 18A calculates Expression 8. 

Y(ni,n2,t) = FFT(ni,n2,X(:,:t)) Expression 8 

In Expression 8, t is the sampling window time index; X(:,:,t) is the complex antenna 
25 matrix (with the colon representing all antenna indexes for one sampling window 
time index); ni is the direction of arrival frequency index; and n2 is the doppler shift 
index. Each FFT calculation is thus a two dimensional FFT calculation on the 
baseband signal, corresponding both to a specific direction of arrival (as depicted by the 
frequency index nO and to a specific doppler shift (as depicted by the frequency index 
30 na). 
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[000183] The output of correlator 50- 18 A, i.e., the Y(ni,n2,t) values computed 
using Expression 8, are stored as correlator output values. The correlator output values 
can be stored, for example, in the correlator output value memory 52-18A of Fig. 18 A. 

[000184] The correlator output analyzer 62-18A of channel estimate (CE) generator 
5 60-18A searches the correlator output values Y(ni,n2,t) and (as step 20-2) determines 
therefrom a maximum absolute value |Y(ni,n2,t)|inax- This maximum absolute value 
|Y(ni,n2,t)|max is utilized by correlator output analyzer 62-18A to determine both the 
direction of arrival (DOA) and time of arrival (TO A) for an arriving wavefront seen in 
the sampling window. In particular, as step 20-3 the correlator output analyzer 62-18A 

10 chooses a sampling window time index t_max at which |Y(ni,n2,t)|niax occurs to be the 
time of arrival of the arriving wavefront. In addition, as step 20-4 the correlator output 
analyzer 62-18A chooses the frequency index ni_max at which |Y(ni,n2,t)|niax occurs to 
determine the direction of arrival (DOA) of the arriving wavefront. Further, as step 20- 
5 the correlator output analyzer 62-18A chooses the index n2_max at which 

15 |Y(ni,n2,t)|max occurs to determine the doppler shift of the arriving wavefront. An 

amplitude for the arriving wavefront is determined as the correlator output analyzer 62- 
18A divides |Y(ni,n2,t)|max by the product of the number of antennas comprising the 
antenna array and the number of sets of pilot data included in the matrix 130 (as step 
20-6). 

20 [000185] Expression 8 and the steps of Fig. 20 represent a generic FFT calculation. 
In a CDMA-specific situation which utilizes a coding generator (such as coding 
generator 30 of Fig. 1), a comparable FFT calculation can be made using a refinement 
of Expression 8 which appears as Expression 9. 

Y(ni,n2,t) = 2 Cj *FFT(ni,n2,X(:,:,t)), j = 1,K Expression 9 

25 Expression 9 is understood from Expression 1, it being further mentioned that Cj is a 
coding sequence symbol value j; and K is a length of the coding sequence. 

[000186] As a result of operation of spatio-temporal joint searcher and channel 
estimator 24- 18 A, an accurate channel estimate can be provided to the detector as a 
spatio-temporal spatial and temporal signature. The spatial signature includes the 
30 direction of arrival; the temporal signature includes the doppler shift. The channel 
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coefficient (CC) for each time of arrival and antenna element is derived from the 
direction of arrival (DOA)and the doppler shift. The time of arrival (TO A) and channel 
coefficient (CC) are applied to the detector as represented by lines 66-18A and 68- 18 A, 
respectively, in Fig. 18A. 

5 [000187] As mentioned above, the channel coefficient (CC) for each wavefront is 
derived from the direction of arrival (DO A) and doppler shift (DS). Recall that at step 
18-4 analyzer 62-18A chose the frequency index ni_max at which |Y(ni,n2,t)|niax occurs 
to represent the direction of arrival (DO A) of the arriving wavefront, with the chosen 
frequency index corresponding to a direction of arrival (e.g., 0). Further, analyzer 62- 

10 18A chose the frequency index n2_max at which |Y(ni,n2,t)|niax occurs to represent the 
doppler shift of the arriving wavefront, with the chosen frequency index corresponding 
to a doppler shift. The channel impulse response vector (i.e., array propagation vector) 
X is therefore generated by detector interface 64-18A in accordance with Expression 10 
(for identical isotropic antenna elements). 

15 X =[(l,e^''^*^^^^^ g,akd*2sine) ^akd*(K-l)sine)-j^:^-.Q. 
Q ^akd*sine) ^akd*2sine) ^Okd*(K-l)sinG)j^^ j . 

(l,e^^^*-^>, e^^^*2sine)^ ^cikd*(K-i)sine)pcN Expression 10 

In Expression 10, CN = d " , with H and other parameters being as previously 
defined. 



20 [000188] In the foregoing description it is the role of channel estimate (CE) 

generator 60- 18 A, and particularly detector interface 64- 18 A, to generate both a time of 
arrival (TO A) and a channel coefficient (CC), the channel coefficient being derived 
from the direction of arrival and doppler shift, e.g., as above described in conjunction 
with Expression 11. In an altemate implementation of this and other embodiments 

25 described herein, the detector itself (such as detector 26 illustrated in Fig. 1), upon 
receiving the time of arrival (TOA), the direction of arrival (DOA), and the doppler 
shift for each arriving wavefront, may have the intelligence to compute the channel 
coefficient for each wavefront from the corresponding direction of arrival (DOA) and 
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doppler shift information. In such case, the time of arrival, direction of arrival, and 
doppler shift are output by detector interface 64 to the detector. 

[000189] The operation of the correlator 50-18A in calculating Expression 8 or 
Expression 9 is an example of a three dimensional essentially concurrent mode, since 
5 evaluation of Expression 8 (or Expression 9 for a WCDM A implementation) involves a 
single step determination of the time of arrival and channel coefficient by simultaneous 
considering signals from all antennas of the array for all of the plural series. In other 
words, in the illustrated example of the three dimensionally essentially concurrent 
mode, the Fast Fourier Transform (FFT) of Expression 8 or Expression 9 had three 
10 arguments: Ui, n2, and X(:,:t), so that the FFT operated on all arguments essentially 
simultaneously. 

[000190] In contrast to the three dimensional essentially concurrent mode, the 
sequenced mode involves a two step determination of the time of arrival and channel 
coefficient. In a first alternative way for implementing the sequenced mode, a first step 

15 comprises determining a time of arrival and direction of arrival by concurrently 

considering the plural signals provided by the plural antennas for a first of the plural 
series. For example, the first step of the first alternative of the sequenced mode can 
involve calculating a FFT such as that of Expression 1 (or, for WCDMA, Expression 
2). From the results of the first step or first FFT calculation, a time of arrival (TO A) 

20 and tentative channel coefficient are determined. Then, as a second step of the first 

alternative of the sequenced mode, the tentative channel coefficient is refined by taking 
into consideration a possible frequency shift (e.g., doppler shift) by further considering 
the elements of the plural series having the direction of arrival determined in the first 
step. In a second alternative way of implementing the sequenced mode, the order of the 

25 steps is essentially reversed: first the FFT is performed in the temporal domain to 
decide time of arrival and tentative channel coefficient; and secondly the tentative 
channel coefficient is refined by FFT in the spatial domain. 

[000191] The procedures of the first alternative implementation of the sequenced 
mode for the non-parametric technique are illustrated in Fig. 22A and Fig. 22B in 
30 conjunction with Fig. 23. Fig. 22A and Fig. 22B as a diagrammatic view of an antenna 
signal matrix; an antenna weight vector; and a non-parametric estimation output vector 
for an example embodiment of a sequential spatio-temporal joint searcher and channel 
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estimator. In Fig. 22A, the FFT operates on the spatial domain and calculates the FFT 
(illustrated by the FFT vector Wi) for the antenna matrix for each time interval. The 
time of arrival is choosen by picking the direction of arrival index and time index with 
the highest absolute value. If this index does not coincide for all time intervals, the 
5 index can be chosen with some method, e.g. majority decision. 

[000192] After having chosen the time of arrival index and direction of arrival 
index, these FFT-processed samples are further FFT-processed by a FFT calculation in 
the temporal domain (illustrated by FFT frequency vector Wj). Fig. 22B shows the 
spatially filtered samples for the identified time of arrival and direction (marked as grey 
10 in the figure) are filtered with the temporal vectors. After the second FFT processing, 
the channel estimate is created from the sample with the highest magnitude. Step 23-1 
through step 23-7 of Fig. 23 also describe the procedure of the first alternative 
implementation of the sequenced mode. 

[000193] The procedures of the second alternative implementation of the sequenced 
15 mode for the non-parametric technique are illustrated in Fig. 24A and Fig. 24B in 
conjunction with Fig. 25. Fig. 24 A and Fig. 24 B show an antenna signal matrix; a 
Doppler weight vector; and a non-parametric estimation output vector. In Fig. 24A the 
FFT operates on the temporal domain and calculates the FFT (illustrated by the FFT 
vector Wj) for the antenna matrix for each time interval. The time of arrival is chosen 
20 by picking the doppler index of and time index with the highest absolute value. If this 
index does not coincide for all time intervals, the index can be chosen with some 
method, e.g. majority decision. After having chosen the time of arrival index and 
doppler index, these FFT-processed samples are further FFT-processed by a FFT 
calculation in the spatial domain (illustrated by FFT frequency vector Wi). Fig. 24B 
25 shows the spatially filtered samples for the identified time of arrival and Doppler shift 
(marked as grey in the figure) are filtered with the spatial vectors. After the second 
FFT processing, the channel estimate is created from the sample with the highest 
magnitude. Step 25-1 through step 25-7 of Fig. 25 also describe the procedure of the 
second alternative implementation of the sequenced mode. 

30 [000194] Whereas the joint searcher and channel estimator of Fig. 1 8 A includes a non- 
parametric type correlator (e.g., a filter which performs a Fast Fourier Transform (FFT) 
calculation), in other example embodiments the joint searcher and channel estimator 
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implements parametric techniques. As does the Fig. 18A embodiment, the parametric 
temporal joint searcher and channel estimator 24-1 SB of Fig. 18B is shown along with 
its associated example antenna array 22-18B. Again by way of example, antenna array 
22-18B includes four antenna elements 22-18B-1 through 22-18B-4. The signals 
5 obtained from the antenna elements are each applied to joint searcher and channel 
estimator 24-18B, as well as to a detector (not illustrated in Fig. 18B). 

[000195] Similarly to the earlier described embodiment, joint searcher and channel 
estimator 24-18B can comprise an antenna signal matrix handling unit 40-18B, which 
in tum comprises antenna signal matrix generator 42-18B and antenna signal matrix 
10 memory 44-1 8B, which function much in the manner previously described. For 

example, the complex baseband values stored in antenna signal matrix memory 44-18B 
can also be conceptualized as matrix 130, and as such has a sampling window time 
index. The antenna signal matrix 80 has been previously discussed in conjunction with 
Fig. 19. 

15 [000196] The joint searcher and channel estimator 24-1 8B further comprises a 

parametric estimator 51-18B which produces a parametric estimation output vector. In 
addition, in similar manner as the preceding embodiment, joint searcher and channel 
estimator 24-18B comprises a channel estimate generator 60-18B which has parametric 
output estimation vector analyzer 62-18B and a demodulator interface 64-18B 

20 [000197] Fig. 26 shows an antenna signal matrix and a parametric estimation output, 
vector for an example embodiment of a spatio-temporal joint searcher and channel 
estimator. As with the non-parametric techniques, the parametric techniques can be 
implemented either in a three dimensional essentially concurrent mode or in a 
sequenced mode, with the sequence mode having two alternative implementations. 

25 [000198] Fig. 27 shows basic, representative steps involved in a parametric three 
dimensional essentially concurrent mode. Step 27-1 shows the joint searcher and 
channel estimator 24-18B producing a parametric estimation output vector. Then, as 
step 27-2, analyzer 62-18B finds the "qualifying" values in the parametric estimation 
output vector. 

30 [000199] For each qualifying value, as step 27-3 the parametric output estimation 
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vector analyzer 62-18B chooses a time of arrival (TOA) as corresponding to the 
sampling window time index t for the qualifying value, e.g., the time index at which the 
maximum/qualifying absolute value of the parametric estimation output vector occurs. 

[000200] For each qualifying value, as step 27-4 the parametric output estimation 
5 vector analyzer 62-18B chooses a spatio-temporal frequency parameter corresponding 
to the spatio-temporal frequency for the maximum/quaUfying absolute value of the 
parametric estimation output vector. 

[000201] As step 27-5, the parametric estimation output vector analyzer 62-13B 
determines the amplitude as the value for the spatio-temporal amplitude value for the 
10 time of arrival decided in step 27-2. 

[000202] It should be understood from the foregoing that information indicative of 
more than one incident wavefront may be seen in a sampling window. For example, 
with reference to the parametric estimate output vector 140 of Fig. 26, the parametric 
output estimation vector analyzer 62-18B may see other (e.g., plural) high numbers, and 
15 for each of those high numbers which qualify, an arriving wavefront may be 
ascertained. 

[000203] The procedures of the first altemative implementation of the sequenced 
mode for the parametric technique are illustrated in Fig. 28A and Fig. 28B in 
conjunction with Fig. 29. Fig. 28A and Fig. 28B depict a parametric, sequential spatio- 

20 temporal joint searcher and channel estimator for this first altemative implementation. 
In Fig. 28A and Fig. 28B the parametric approach first operates on the spatial domain 
and calculates the spatial frequency parameters for each time instant over the time 
transmission intervals. The time of arrival is chosen by picking the spatial frequency 
amplitude value with the highest absolute value. The direction of arrival, DOA, is the 

25 value of the spatial frequency parameter. If this time of arrival does not coincide for all 
the time intervals, the time of arrival can be chosen with some method, e.g. majority 
decision. As shown in Fig. 28B, after having chosen the time of arrival index and 
direction of arrival, these samples are processed by the parametric approach applied in 
the temporal domain. After the second processing, the channel estimate is created from 

30 the temporal parameters. Step 29-1 through step 29-5 of Fig. 29 also describe the 
procedure of the first altemative implementation of the parametric sequenced mode. 
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[000204] Fig. 30A and Fig. 30B show a parametric, sequential spatio-temporal 
joint searcher and channel estimator for a second alternative implementation of a 
parametric, sequential spatio-temporal joint searcher and channel estimator. In Fig. 
30A and Fig. BOB the paremetric approach first operates on the temporal domain and 
5 calculates the temporal frequency parameters for each time instant over the time 

transmission intervals. The time of arrival is chosen by picking the temporal frequency 
amplitude value with the highest absolute value. The Doppler shift frequency, DSF is 
the value of the temporal frequency parameter. If this time of arrival does not coincide 
for all the time intervals, the time of arrival can be chosen with some method, e.g. 
10 majority decision. As shown in Fig. BOB, after have chosen the time of arrival index 
and DSF, these samples are processed by parametric approach applied in the spatial 
domain. After the second processing, the channel estimate is created from the spatial 
parameters. Step Bl-1 through step Bl-7 of Fig. Bl also describe the procedure of the 
second alternative implementation of the parametric sequenced mode. 

15 [000205] The non-parametric FFT-type correlator and the parametric linear 

combination logic techniques illustrated above are only two example techniques for 
finding the values or "tones" in antenna signal matrix 130 which are associated with 
arriving wavefronts. Other parametric approaches are described by or understood from 
Stocia, Petre and Moses, Randolph, Introduction To Spectral Analysis^ ISBN-013- 

20 258419-0, Prentice Hall, which is incorporated by reference in its entirety, particularly 
Chapter 4 thereof. 

[000206] The spatio-temporal joint searcher and channel estimator and techniques 
of operation thereof as described above are suitable for any receiver unit which has 
plural receiving antennas. Thus, the spatial joint searcher and channel estimator is 
25 particularly well suited for, but not limited to, a base station which has plural antennas. 
The spatio-temporal joint searcher and channel estimator and techniques of operation 
thereof also encompasses mobile terminals which have plural antennas. 

*** 

[000207] The joint searcher £ind channel estimators thus employ a multi- 
30 dimensional and optimum detection and estimation approach. The multi-dimensional 
joint searcher and channel estimators typified by those described herein have better 
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performance than traditional one dimensional searchers. The multi-dimensional joint 
searcher and channel estimators have a greater SNIR for detecting time of arrival, 
which increases the probability of that the correct time of arrival will be ascertained. 
This, in turn, leads to a better channel estimate. 

5 [000208] In terms of implementation, the blocks, units, and functionalities of the 
differing embodiments of the joint searcher and channel estimator herein described can 
various forms. For example, those skilled in the art will appreciate that one or more of 
the functionalities of the joint searcher and channel estimator can be implemented using 
individual hardware circuits, using software functioning in conjunction with a suitably 

10 programmed digital microprocessor or general purpose computer, using an application 
specific integrated circuit (ASIC), and/or using one or more digital signal processors 
(DSPs). Furthermore, the functionalities of the joint searcher and channel estimator 
need not be delineated specifically in the manners illustrated, it being understood (for 
example) that the functionaUties can be distributed, combined, subdivided, or otherwise 

15 rearranged for accomplishing essentially the same results. 

[000209] Use and operation of the joint searcher and channel estimators is not 
confined to WCDMA transmission, although in some instances WCDMA has been 
described above as an example environment of implementation. The principles, 
techniques, methods, and apparatus described herein can be adapted or augmented for 
20 compatibility with various types of networks, not only WCDMA, but other networks as 
well (such as GSM, for example). 

[000210] In the foregoing, it will be appreciated that other aspects of wireless 
receiver structure and operation which are tangential to matters described above have 
been omitted for clarity. Such aspects, well understood by persons skilled in the art, 
25 include without limitation pulse shaping, sampling frequency, time jitter, time 

alignment, demodulation, inter symbol interference (ISI), and co channel interference 
(CCI). 
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[00021 1] While the invention has been described in connection with what is 
presently considered to be the most practical and preferred embodiment, it is to be 
understood that the invention is not to be limited to the disclosed embodiment, but on 
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the contrary, is intended to cover various modifications and equivalent arrangements 
included within the spirit and scope of the appended claims. 



